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The  general  objective  of  the  studies  conducted  by  the  Agricultural  Research  Service,  Ithaca,  N.Y., 
is  to  evaluate  plant,  soil,  and  meteorological  interactions  involved  in  the  partition  of  thermal  energy  at 
the  earth’s  surface.  This  report  presents  in  part  the  results  of  research  conducted  with  the  particular 
objective  of  evaluating  the  aerodynamic  surface  roughness  and  relating  it  to  the  elastic  and  geometric 
characteristics  of  the  surface  cover. 

Theoretical  and  experimental  investigations  of  the  turbulent  iransfer  characteristics  of  the  air- 
stream  near  the  ground  are  reported.  It  is  shown  that  the  surface  boundary  layer  must  he  divided  into 
two  regions:  the  freestream  above  the  surface  where  the  various  forms  of  turbulent  transfer  are  nearly 
independent  of  height,  and  the  airstream  within  the  vegetative  canopy  where  sources  and  sinks  are 
present.  Several  theoretical  models  for  canopy  flow  are  reported  ami  compared  with  ohi  ovations. 
A  numerical  method  for  determining  the  aerodynamic  characteristics  and  their  standard  errors  in  the 
freestream  is  described.  Experimental  results  for  freestream  and.  canopy  flow  are  reported. 
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model  in  which  the  mixing  length  is  a  function  of  height  above  the  ground  surface.  Experimental  data 
show  that  the  aerodynamic  surface  roughness  varies  with  the  windspeed  and  the  variations  depend  upon 
the  geometric  and  elastic  characteristics  of  the  vegetation.  Elastic  vegetation  shows  decreasing  rough¬ 
ness  with  increasing  windspeed  and  shearing  stress  is  nearly  independent  of  windspeed.  Semirigid 
vegetation  has  increasing  roughness  and  shearing  .dress  with  increasing  windspeed.  Tin  results  of  these 
studies  serve  to  emphasize  that  the  concept  of  “roughness  length,”  although  very  useful,  is  in  reality  an 
artifact. 
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THE  ENERGY  BUDGET  AT  THE  EARTH’S  SURFACE 

Part  H 

INTRODUCTION 


Unfortunately,  the  evaluation  of  the  partition 
of  solar  radiation  into  other  forms  of  energy  at  the 
earth’s  surface  has  to  depend,  in  part,  upon  indirect 
means.  For  example,  energy  exchange  calcula¬ 
tions  based  upon  aerodynamic  principles  have  to 
be  used  to  estimate  sensible  heat  exchange.  The 
principles  are  customarily  used  also  to  estimate 
latent  heat  and  momentum  exchange.  More  re¬ 
cently,  the  photochemical  fixation  of  carbon  di¬ 
oxide  in  the  process  of  photosynthesis  has  been 
evaluated  from  aerodynamic  exchange  calculations 
of  carbon  dioxide  transfer  to  vegetated  surfaces. 

Historically,  the  aerodynamic  exchange  princi¬ 
ples  were  first  developed  in  wind  tunnel  and  lab¬ 
oratory  fluid  flow  studies,  then  applied  to 
atmospheric  exchange  studies  at  the  earth’s  surface 
either  over  water  or  over  land  where  the  vegetation 
was  generally  short  or  sparse.  It  has  been  com¬ 
monly  accepted  that  such  calculations  are  more 
complicated  over  water  because  of  its  surface  in- 
stabdity.  This  arises  from  the  fact  that  in  the 
exchange  calculations  commonly  used,  it  is  neces¬ 
sary  to  know  certain  boundary  conditions,  i.e., 
have  knowledge  of  a  reference  plane  where  the 
turbulent  boundary  layer  begins.  One  can  readily 
see  that  on  land  surfaces  the  problem  should  be 
much  simpler,  since  the  stable  soil  surface  can  serve 
as  the  reference  plane. 

Recent  investigations  over  more  lush  vegetative 
surfaces  of  the  land,  however,  have  demonstrated 
that  these  surfaces  may  theoretically  and  experi¬ 
mentally  be  more  difficult  to  handle  than  water 
surfaces  simply  because  of  changing  geometric 
characteristics  (plant  growth,  leaf  fall)  and  chang¬ 
ing  elastic  properties  (stalk  bending,  leaf  flutter, 
and  streamlining;  Of  course,  these  problems  are 
amplified  with  larger  plants;  i.e.,  prairie  grass, 
farm  crops,  and  forests. 

When  one  wants  to  evaluate  the  heat  budget 
over  extended  periods  of  time— for  days  and 
months  as  we  are  attempting  to  do — it  is  prac¬ 
tically  impossible  at  present  to  use  any  aerody¬ 
namic  technique  other  than  the  so-called  log 
rofile  methods.  Since  this  method  depends  upon 
nowledge  of  the  surface  reference  plane  where 
the  turbulent  boundary  layer  begins  (and  we  are 
particularly  interested  in  the  energy  exchanges 


with  extensive  vegetative  surfaces),  it  is  necessary 
to  resolve  the  new  complications  that  have  nat¬ 
urally  developed.  This  requires  a  systematic 
research  effort  into  the  turbulent  exchange  charac¬ 
teristics  within  and  above  the  dynamic  vegetative 
canopy.  This  report  takes  up  in  detail  the  theo¬ 
retical  and  experimental  approaches  used  in  1960 
to  understand  further  the  aerodynamic  exchange 
processes  within  and  above  the  plant  canopy. 
The  studies  were  carried  out  in  addition  to  t&e 
routine  collection  of  heat  budget  data  for  the  1960 
growing  season  over  an  alfalfa  field  and  a  cornfield. 
The  results  of  the  heat  budget  studies  will  be 
taken  up  in  another  publication. 

The  nrst  three  sections  in  the  present  study 
take  up  the  theoretical  development  and  support¬ 
ing  experimental  results  concerning  the  airflow 
characteristics  within  the  source  and  sink  region 
of  the  turbulent  boundaiy  layer  at  the  earth’s 
surface,  as  well  as  the  airflow  characteristics  just 
above  the  source  an.,  dink  region.  Two  models 
for  ‘‘canopy  flow”  (in  the  source  and  sink  layer) 
were  developed.  It  was  later  found  that  the 
model  based  upon  fully  turbulent  flow  within  the 
plant  canopy  most  realistically  applies  to  nature. 
To  facilitate  experimental  studies,  it  was  postu¬ 
lated  that  the  turbulent  structure  could  be  broken 
down  into  two  classifications:  (1)  a  transient  state 
and  (2)  a  quasi-steady  state.  To  date,  only  the 
quasi-steady  state  has  been  investigated.  The 
turbulent  canopy  flow  model  postulated  a  linear 
growth  of  diffusivity  within  the  canopy  flow  layer, 
rather  than  the  log  velocity  profile  (which  holds 
above  the  canopy  flow  layer).  The  theory  and  ex¬ 
perimental  results  demonstrate  a  curvature  re¬ 
versal  of  the  velocity  profile  deep  in  the  canopy 
layer. 

The  last  two  sections  discuss  the  complexity  of 
the  relationship  of  the  logarithmic  velocitv  profile 
characteristics  above  the  canopy  to  the  dynamic 
surface  properties  of  lush  vegetation.  To  facili¬ 
tate  the  study  of  the  relationships  involved,  a 
machine  computation  method  was  worked  out  and 
described  (p.2K).  From  the  method  one  can  eval¬ 
uate  the  friction  velocity,  roughness  length,  and 
effective  displacement  parameters  along  with  their 
standard  errors.  The  last  section  takes  up  the 
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experimental  findings  from  log  profile  studies  over 
alfalfa,  wheat,  and  corn.  An  amazing  phenom¬ 
enon  was  discovered  in  that  both  the  alfalfa  and 
wheat  demonstrated  a  constant  friction  velocity 
with  increasing  windspeed  once  a  “critical”  win<i- 
speed  was  reached.  Evidently,  the  elastic  plants 
“streamlined”  so  that  a  decrease  in  roughness 
compensated  for  increasing  windspeed. 

The  last  section  also  takes  up  an  evaluation  of 


the  distribution  of  shear  and  the  momentum  trans¬ 
fer  coefficient  within  the  vegetative  canopy  in 
wheat  and  corn.  These  studies  focus  attention 
on  the  fact  that  shear  originates  within  the  rough¬ 
ness  volume,  and  not  at  some  empirical  plane. 
Further  studies  are  to  be  made  to  evaluate  also 
the  source  and  sink  distribution  for  sensible  heat, 
latent  heat,  and  carbon  dioxide  within  the  canopy 
layer. 


EFFECTS  OF  TURBULENCE  AND  PHOTOSYNTHESIS  ON  CO 
PROFILES  IN  THE  LOWER  ATMOSPHERE 

By  D.  E.  Ordway,  A.  Ritter,  D.  A.  Spence,  and  H.  S.  Tan 


BACKGROUND 

The  two  causes  of  turbulence  in  the  atmosphere 
are  (1)  forced  convection,  through  the  boundary 
layer  type  shear  flow  instability  (effect  of  velocity 
profile),  and  (2)  free  convection,  through  the  buoy¬ 
ancy  force-type  instability  (effect  of  heating  from 
below). 

At  present  much  literature  on  the  study  of 
turbulence  exists  on  forced  convection,  but  little 
definitive  knowledge  is  available  regarding  the 
second  type.  In  the  present  study,  it  is  our  pur¬ 
pose  to  limit  the  first  phase  of  the  investigation 
to  the  effect  of  the  forced  convective,  or  boundary- 
layer  type  turbulence  only,  on  atmospheric  Cl)2 
distribution  We  are  considering  the  problem  of 
steady  state  atmospheric  boundary  layer  flow  and 
shall  now  introduce  a  principle  of  “flux  conserva¬ 
tion.” 

This  principle  states  in  essence  that  the  trans- 

Eort  equation  for  any  physical  quantity  Q  must 
e  of  the  form 

^=^+g.V£=div(iev£)+SsO  (l) 

where  K  is  the  flux  coefficient,  S  the  distributed 
source  distribution,  t  is  time,  and  q  the  flow 
velocity.  In  other  words,  particle  derivative  of 
Q  must  vanish.  This  is  indeed  a  generalization  of 
the  well-known  constant  flux  hypothesis  (29)  1 
long  adopted  by  micrometeorologists,  to  include 
the  distributed  source  contribution. 

Under  the  usual  hypothesis  of  flow  stratifica¬ 
tion,  application  of  this  principle  to  the  steadv 
state  distribution  of  C()2  concentration,  E,  witli 
height,  z,  immediately  gives 

v.  pV/.;=|z[Dd|]+«,.0  (2) 

where  ir  is  distributed  sink  or  source  of  (X)2  due 
to  photosynthesis  and  respiration,  and  l)  the  total 
diffusivity;  i.e.,  suin  of  molecular  diffusion  (lam¬ 
inar,  Di)  and  eddy  mixing  (turbulent,  l)r). 

Application  of  the  same  principle  to  x-moiiien- 
tum  with  w  as  wind  velocity  gives 

q.  PV»=^-d  =  ()  (3) 

uz 


where  d  is  distributed  shear  source,  i.e.,  leaf  drag, 
and  r  is  the  shear  stress,  given  by 


with  p  denoting  laminar  viscosity  and  «,  eddy 
viscosity. 

The  validity  of  this  principle  can  most  easily 
be  seen  from  the  limiting  behavior  of  equation  3. 
Indeed,  for  outer  flow  where  rfsO,  this  equation 
leads  to  a  constant  shear  requirement  that  is 
famous  for  the  logarithmic  velocity  profile  in 
turbulent  wall  flow  and  linear  velocity  profile  in 
laminar  f'ouette  flow.  On  the  other  hand,  when 
this  equation  directly  gives  the  differential 
equation  of  canopy  flow  velocity  distribution. 

Equations  2  and  3  will  form  the  basis  of  our 
mathematical  study  of  the  atmospheric  carbon 
dioxide  concentration  and  the  canopy  flow,  respec¬ 
tively. 

ATMOSPHERIC  CARBON  DIOXIDE 
CONTENT  STUDY 

COj  content  in  the  lower  atmosphere  controls 
and  is  in  turn  controlled  by  the  life  process  on 
earth.  The  mechanism  that  keeps  local  C02  con¬ 
tent  from  becoming  deficient  or  excessive  is  the 
process  of  molecular  (laminar)  and  eddy  (turbu¬ 
lent)  diffusion.  The  differential  equation  gov¬ 
erning  the  steady  state  COj  concentration,  E,  as¬ 
suming  adequate  (X)2  for  saturation  requirement 
of  photosynthesis,  can  therefore  be  written  in 
general  as  follows: 

bz  (D1Tz)+W=0  (5) 

w=—j(z)E+g(z)  (6) 

where  V  denotes  the  total  diffusivity,  w  the  rate  of 
('02  assimilation  or  production  as  a  result  of  photo¬ 
synthesis  (/)  and  respiration  ( g ).  Functions  j  and 
g  are  defined  by: 

f(z)  =  photosvnthetie  rate  coefficient^ \L(z)c'(z) 

_  '  (7) 

1  Italic  figures  in  parentheses  refer  to  Literature  Cited, 
p.  -ts. 
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<7(2)— rate  of  respiration =Cj»' (2).  (8) 

Here  L(z)  is  the  light  intensity  function  (5),  a'(z) 
the  leaf  area  density  function,  and  C\  and  C2  are 
constants. 

With  tt’=0  for  2>A  (^=plant  height),  it  is 
evident  that  the  plane  z=~h  separates  the  atmos- 

f there  into  two  regions,  inside  each  of  which  dif- 
erent  laws  govern  the  C02  distribution,  namely. 

The  outer  flow,  z^>h — 

(9> 

The  inner  flow,  2<A — 

Up\ !)-/WE+?<d=0-  (10) 

Equations  5  and  6  show  clearly  the  respective 
roles  played  by  turbulence  and  photosynthesis  on 
the  atmospheric  C02  distribution:  the  effect 
of  turbulence  appears  in  the  diffusive  term 

through  D,  the  eddy  diffusivity  coeffi¬ 
cient.  The  effect  of  photosynthesis  that  depends 
on  leaf  area,  light  intensity,  and  atmospheric  C02 
concentration  appears  as  the  linear  term  j{z)E  of 
the  differential  equation.  The  nonhomogeneous 
term  g(z)  represents  the  effect  of  respiration,  which 
is  independent  of  light  and  C02  contents. 

Now,  the  kinetic  transport  mechanism  for  mass 
(diffusion),  momentum  (shear),  and  energy  (heat 
conduction)  are  identical.  Upon  introducing  the 
concept  of  constant  turbulent  Schmidt  number 
defined  by 


01) 


we  can  usually  establish  some  correlation  between 
C02  concentration  profile  E(z)  and  velocity 
profile  u{z).  Let  us  consider  now  the  two  flow 
regimes. 


However,  one  must  be  careful  in  not  pushing 
the  analogy  too  far.  Equation  12  is  indeed 
identical  with  that  governing  the  velocity  dis¬ 
tribution,  but  the  analogy  stops  here.  Wnereas 
the  nonslip  condition  over  the  datum  plane  of 
velocity  profile  necessarily  represents  a  momentum 
sink,  the  boundaiy  condition  over  the  datum 
plane  of  concentration  profile  can  be  either  source 
or  sink,  depending  solely  on  the  relative  rate  of 
respiration  and  photosynthesis.  Moreover,  even 
the  datum  planes  for  the  two  cases  are  usually  at 
different  locations.  Thus,  the  only  conclusion 
one  can  safely  draw  from  the  similitude  of  the  two 
differential  equations  is  that  the  two  solutions 
must  belong  to  the  same  family  of  curves. 

Inner  Flow,  z</» 

In  this  region,  r  ^constant.  As  can  be  seen  by 
comparing  with  the  canopy  flow  equation  derived 
in  the  discussion  of  the  laminar  flow  model  study, 
p.  5,  the  differential  equation  for  E  differs 
essentially  from  that  for  u.  Thus  little  correlation 
between  the  two  profiles  can  be  expected.  Here, 
two  cases  should  be  considered  separately,  de¬ 
pending  on  the  linear  term  in  the  differential 
equation. 

1.  Respiration  Profile  ( Nighttime ). 

In  this  case, 

/(z)=Z(z)-O,0(z)*O.  (14) 

Equation  10  reduces  to  the  following  form, 

|(cg)+^)=°.  (1« 


The  solution  to  this  equation  can  be  obtained 
through  integration,  as  follows: 


E— 


dE 

dz 


=5  [^- J</(2)^] 
E0=j  7)  |/1— 


(16) 


Outer  Flow,  z>A 

Since  outer  flow  is  nothing  but  the  turbulent 
wall  flow,  where  t  =  t0,  t  —  Cz,  on  introducing  the 
Schmidt  number,  equation  9  takes  the  form 


giving 

E-E0=C\nz.  (13) 

Like  the  velocity  distribution,  E  also  has  a 
logarithmic  profile. 


2.  Photosynthesis  Profile  {Daytime). 

In  this  case,  if 


.f(z)* 0,  g(z)y* 0. 

the  full  equation  10  must  be  retained;  i.e., 


D 


d*E  dD  dE 
dz 1  dr  dz 


-f(z)E——g(z). 


(17) 

(18) 


A  general  solution  to  this  equation  usually  pre¬ 
sents  great  complication  and  difficulty.  How¬ 
ever,  if  we  make  the  approximation  D=constant, 
which  has  been  commonly  adopted  by  various 
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where  n  is  the  laminar  viscosity  and  F(z)  is  some 
function  of  the  plant  characteristics.  The  indi¬ 
rect  solution  to  this  equation  can  be  obtained  by 
quadrature  after  certain  transformations  and  ex¬ 
pressed  in  terms  of  elliptic  functions  as  follows: 


Although  this  assumption  is  admittedly  crude, 
nevertheless,  it  should  provide  a  good  first  ap¬ 
proximation. 

A  solution  to  equation  19  can  be  obtained  by 
introducing  a  Fourier  spectrum  for  j(z)  and  ex¬ 
pressing  the  complementary  solution  in  the  form 


F^u^u2du+^ 

f  “  du 

2  Jo  lQ(u)+kl)',2+  2' 


investigators  (32),  this  equation  simplifies  to  the 
following  generalized  Hille’s  equation: 

Dj?-j(z)E=-g(z).  (19) 


Y=ea^ane'D‘+ea2‘^2bne'“=Yi  +  Y2.  (20) 

The  determination  of  characteristic  exponents  a, 
by  Hille  determinant  and  the  recursion  formulas 
connecting  the  an’ s  and  6,’s  are  discussed  in  detail 
by  Whittaker  and  Watson  (39).  For  the  partic¬ 
ular  solution,  by  writing 

Y—V\Yi=v2Y2  (21) 

and  applying  the  method  of  variation  of  para¬ 
meters,  it  can  be  shown  that  (10) 

Vl=~fl(YuY7)  9{2)d2 

t'2=/A(rI,F2)  9^z)dz  (22) 

where  A(FhF2)  is  the  Wronskian  of  F]  and  F2. 


A  detailed  discussion  and  some  typical  results  of 
this  model  are  given  below. 

In  the  first  model  of  the  canopy  flow  velocity 
distributions  study,  an  approximate  solution  to 
the  velocity  profile  u(z)  in  tne  base  flow  within  the 
plants  has  been  considered.  At  the  upper  edge 
of  the  flow,  say  z  =  z ,,  we  see  that  the*  outside 
turbulent  flow  effectively  “drags”  the  canopy 
flow,  analogous  to  a  moving  wall  with  a  velocity 
u,.  If  we  replace  the  leaves  by  an  effective  distri¬ 
bution  of  infinitesimal  “shear  centers,”  we  can 
write  an  equation  of  horizontal  force  equilibrium 
from  figure  1.  We  assume  that  r,  the  fluid  dy¬ 
namic  shear,  is  due  mainly  to  the  molecular 
momentum  transport  of  laminar  flow  and  D,  the 
average  drag  per  unit  volume  of  the  shear  centers, 
is  turbulent  in  nature.  In  nondimensional  form, 
this  yields, 

0HS*“S(7)  f  (25) 


CANOPY  FLOW  VELOCITY 
DISTRIBUTION  STUDY 

We  define  canopy  flow  as  the  flow  inside  the 
plants  (0 <2<A),  which  consists  mainly  of  a 
turbulent  layer  at  the  upper  part  and  a  wake-type 
flow  at  the  lower  part.  It  is  clear  that  a  ready 
satisfactory  study  must  be  a  statistical  one,  of 
ensemble  average  type  approach.  However,  as  a 
crude  approximation,  it  seems  also  plausible  that 
we  can  replace  the  distributed  plant  leaves  by  a 
fictitious  distribution  of  shear  centers  and  solve 
the  differential  equation  corresponding  to  this 
model.  This  approach  is  by  nature  speculative. 
At  the  present  moment  there  are  two  models  of 
approach  under  consideration. 

Laminar  Flow  Model 

The  differential  equation  for  velocity  u  can  be 
written: 

^-F(z)u2=0  (23) 


where  u  *  u/u , ,  z  =•  2/2 ,  and 

S(z)  ssCD  Ret  (26) 

CD  is  a  mean  drag  coefficient  of  the  shear  centers, 
which  is  turbulent  and  is  independent  of  u;  Ret  is 
the  Reynolds  number  dependent  on  ut  and  z,,  and 
d<r(z)laz  is  an  effective,  equivalent  local  leaf  area 
density.  Equation  25,  together  with  the  boundary 
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conditions  at  z,  and  0,  respectively, 
5(1)  =  1 
5(0)  =0 


Further  study  of  both  models  should  be  carried 
out.  A  Anal  choice  can  only  be  made  after 
(27)  sufficient  field  data  have  been  gathered  to  confirm 
the  theoretical  model. 


then  determines  the  velocity  profile  in  the  canopy 
flow.  Since  the  solution  is  characterized  only  by 
S(z),  we  call  S(z)  the  “shape  function.”  Because 
it  is  always  positive,  we  see  that  the  variation  of 
the  slope  is  always  negative.  That  is,  the  velocity 
profile  is  always  concave  backwards  as  opposed 
to  the  typical  laminar  boundary  layer,  which  is 
concaved  forward. 

Equation  25  is  a  auasi-nonlinear  differential 
equation  of  second  order.  In  general,  the  solu¬ 
tion  has  not  been  found.  However,  for  constant 
S  or  SCz)=5,  a  suitable  transformation  to  the 
velocity  gradient  as  the  independent  variable  and 
the  velocity  as  the  dependent  vaiiable  permits 
integration.  The  result  is 


:=f£cos‘ 


(V5—  l)k— u 
(V3+l)*+« 


-f[cos-*(2-V5),  (29) 


where  k  is  related  to  S  and  the  profile  slope  at  the 
ground;  i.e.,  the  ground  shear  by 


*4(1)1 


(30) 


and  F  is  the  elliptic  function  of  the  first  kind. 
For  a  given  value  of  the  constant  shape  function, 
the  parameter  k,  and  hence  the  ground  shear,  is 
determined  in  general  by  satisfying  equation  27. 
The  results  are  plotted  in  figure  2  over  a  repre¬ 
sentative  range.  A  typical  velocity  profile  for 
5=141  has  also  been  computed  and  is  illustrated 
in  figure  3. 

Turbulent  Flow  Model 


This  model  has  been  suggested  under  the  heu¬ 
ristic  argument  that  on  the  leaf  scale,  the  flow  may 
appear  laminar;  on  plant  scale,  it  must  appear 
turbulent.  The  corresponding  differential  equa¬ 
tion  takes  the  form: 


i('d£)-Fw=°  <3» 

where  t(z)  is  the  eddy  viscosity. 

For  a  first  approximation,  we  may  take  «  as 
constant  (several  orders  of  magnitude  greater 
than  n  (S3)),  thus  reducing  equation  31  to  equation 
24,  only  with  a  different  scale. 


k 


ii/4i  )q 


Figure  2. — Canopy  flow  wall  shear  for  a  constant  shape 
function. 


Relative  Haight  (i) 


Figure  3. — Canopy  flow  velocity  profile. 


QUASI-STEADY  MICRO-METEOROLOGICAL  ATMOSPHERIC 
BOUNDARY  LAYER  OVER  A  WHEATFIELD 


By  H.  S.  Tan  and  S.  C.  Ling 


THEORETICAL  MODEL 

Concept  of  Relevant  Scale 

In  the  study  of  turbulent  transport  process  in 
the  lowest  part  of  the  earth’s  atmosphere,  one  of 
the  most  important  concepts  is  tne  “relevant 
scale”  of  the  micro-structure  of  the  wind.  Indeed, 
anv  typical  wind  record  taken  at  a  point  in  a  corn¬ 
field  or  wheatfield  exhibits  (fig.  4)  at  least  two  dis¬ 
tinct  ranges  of  eddy  scales:  (1)  a  high  frequency 
fluctuating  velocity  component  of  small-scale 
eddies  superimposed  on  (2)  a  lower  frequency 
gusty  wind  composed  of  large  eddies.  To  a  cosmic 
observer  in  the  large  or  meteorological  sense,  both 
these  eddies,  with  scales  of  different  orders  of  mag¬ 
nitude,  are  nothing  but  turbulence.  However,  to 
an  observer  sitting  at  the  top  of  a  wheat  crop  in 
the  small  or  micro-meteorological  sense,  the  pic¬ 
ture  is  different.  What  he  recognizes  here  is  a 
turbulent  (high  frequency),  gusty  (low  frequency) 
wind.  The  recognition  of  such  a  "relevant  scale” 
is  important  if  we  recall  that  the  term  “turbulence” 
is  so  vaguely  defined.  Indeed,  P.  A.  Sheppard 
brought  this  point  out  most  clearly  in  his  final 
address  to  the  International  Symposium  on  At¬ 
mospheric  Turbulence  in  the  Boundary  Layer  (7). 
He  points  out  as  follows: 

The  atmosphere  .  .  .  shows  simultaneously  many 
modes  of  motion,  varying  in  scale  from  the  “long 
waves”  of  the  westerlies  down  through  depressions 
and  anticyclones,  convective  motions  and  so  on,  to 
small  scale  turbulence  .  .  . 

Now  the  motions  which  we  commonly  study  under 
the  name  of  turbulence  are  so  much  smaller  in  scale  than 
depressions  and  the  like  that  we  make  little  or  no 
attempt  to  pick  out  individual  elements  of  the  flow, 
except  quite  formally  as:  c.g.,  in  a  power  spectrum. 
But  who  doubts  that  if  we  magnified  the  field  suffi¬ 
ciently,  drew  a  series  of  “synoptic  charts,”  in  fact, 
we  should  in  general  find  recognizable  structures  and 
changes  in  those  structures  with  time  .  .  . 

If  we  study  the  wind  structure  in  its  spectrum 
over  a  relatively  short  time,  we  notice  that  the 
energy  of  turbulence  usually  resolves  itself  into 
two  widely  separated  wave  packets  (fig.  5).  This 
phenomenon  can  be  interpreted  from  a  micro-mete¬ 
orological  point  of  view  as  follows:  Eu  associated 
with  small-scale  eddies,  carries  the  turbulent 
mixing  part  of  the  kinetic  energy  that  controls 
both  the  vertical  and  horizontal  transport  proc¬ 
esses;  Ei  associated  with  large-scale  ecfdies, 


Tim*  (sac.) 


Figure  4. — Horizontal  wind  (u)  at  7-metcr  height. 

IU) - 
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X  (maters) 

Figure  5. — Typical  spectrum  of  horizontal  wind. 


represents  the  part  of  the  kinetic  energy  carried 
away  by  gusty  winds,  which  can  only  affect  the 
horizontal  transport.  This  separation  of  two  dis¬ 
tinct  eddy  scales  responsible  tor  the  vertical  and 
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the  horizontal  transport  processes,  respectively, 
was  also  arrived  at  by  A.  Monin  (20). 

Subdivision  of  Atmospheric  Boundary  Layer: 

Synoptic  Versns  Micro-Layers 

Whenever  there  is  a  fluid  flow  past  a  surface,  a 
boundary  layer  will  develop.  Now  the  structure 
of  the  atmospheric  wind  consists  of  components 
of  large-  and  small-scale  eddies  embedded  in  the 
synoptic  (geostrophic)  wind.  The  large-scale 
turbulence  present  in  the  synoptic  wind  boundary 
layer  is  caused  by  the  presence  of  large-scale 
roughness  elements,  such  as  terrain  irregularities, 
and  the  layer  thickness  usually  extends  to  a  height 
of  a  few  kilometers.  The  relatively  smaller  scale 
micro-meteorological,  or  gusty  winds,  are  nothing 
but  these  large-size  turbulent  eddies  in  the  synoptic 
layer.  Indeed  these  are  actually  air  lumps  of  the 
size  of  a  few  hundred  meters,  sweeping  across  the 
field  and  leaving  behind  the  gust  front  a  growing 
micro-boundary  layer.  This  layer  quickly  reaches 
a  quasi-stndy  stale  a  few  meters  behind  the  prop¬ 
agating  front  and  resembles  very  much  the  steady 
flow  leading  edge  boundary  layer  over  a  flat  plate. 
The  thickness  of  this  micro-boundary  layer  inside 
the  micro-meteorological  wind  usually  extends  to 
about  10  meters.  Further,  the  small-scale  tur¬ 
bulence  inside  this  layer  with  eddy  sizes  ranging 
in  decameters  is  caused  by  small-scale  roughness 
elements,  such  a  plant  crops. 

On  the  average,  a  gusty  wind  of  10  meters  per 
second  with  a  steady  interval  of  several  minutes 
and  a  turbulent  fluctuating  component  having  a 
frequency  of  a  fraction  of  a  cycle  per  second  are 
typical.  On  the  other  hand,  a  micro-wind  in  a 
field  is  never  steady.  However,  the  relatively 
well-defined  edge  of  gusts  (air  lumps)  indicates 
that  a  frontal  theory  can  be  developed.  The 
highly  transient  portion  behind  the  gust  front 
extends  about  10  meters,  after  which  a  quasi¬ 
steady  state  is  quickly  established.  This  means 
that  a  10-second  average  field  measurement  taken 
in  the  middle  portion  of  a  1-ininute  period  gusty 
wind  should  indicate  a  consistent  steady  state 
logarithmic  velocity  profile  outside  the  crop. 
The  validity  of  this  scaling  idea  has  now  been 
amply  verified  by  field  measurements. 

The  classification  of  eddy  sizes  and  definition  of 
a  quasi-steady  condition  can  be  quantitatively 
expressed  in  terms  of  the  micro-boundary  layer 
thickness,  6,  and  the  relaxation  time,  (,  in  the 
following  manner. 

Quasi-steady  condition  considered  reached  at : 
where  6=  10  m., 

t)>Jr  sec. 
c  0 

<=0  at  the  gust  front. 


Now  let  X  denote  the  eddy  size,  or  wavelength, 
and  then  we  have  the  following  scheme: 


Scale 

Eddy  clatti fication 

?<20-- . 

_ Turbulent  eddies. 

2.0<v<20 . 

0 

_ Unsteady  gusts. 

J>20 . 

_ Quasi-steady  wind. 

Proper  Averaging  Period 

In  a  natural  field,  over  a  sufficiently  long  period, 
all  scales  of  eddies  will  appear  and  the  spectrum 
will  be  continuous.  A  velocity  mean  taken  over 
an  arbitrary  period  of  time  is  thus  meaningless. 
The  proper  averaging  period  must  be  arrived  at 
by  “relevant  scale”  conside.  .tions  in  accordance 
with  the  process  involved.  For  example,  in  the 
vertical  transport  process  that  is  controlled  solely 
by  small-scale  turbulence  (wavelength  of  meter 
size  and  period  of  seconds),  extraneous  effects  of 
the  large-scale  turbulence  (gusty  winds)  will  be 
introduced  if  the  averaging  period  is  taken  in 
hours. 

This  will  utterly  ruin  the  statistical  behavior  of 
the  relevant  small-scale  eddies.  On  the  other 
hand,  for  the  horizontal  transport  process  that  is 
governed  by  all  scales  of  turbulence,  averaging 
periods  of  less  than  a  minute  will  give  no  trend  at 
all.  Moreover,  the  results  will  be  meaningless 
since,  under  such  an  averaging  process,  the  largest 
energy-carrying  eddy  scales  in  the  spectrum  are 
completely  left  out. 

The  technique  proposed  in  the  present  report 
therefore  involves  a  selective  sampling  procedure, 
in  order  to  insure  the  proper  quasi-steady  be¬ 
havior.  It  must  be  remarked  at  this  point  that 
this  quasi-steady  state  study  alone  is  not  sufficient 
for  complete  knowledge  of  the  vertical  transport 
process;  this  must  be  augmented  by  a  further  study 
of  the  transient;  i.e. ,  an  accelerating  and  decelerat¬ 
ing  regime.  The  final  effective  eddy  diffusivity 
will  have  to  be  evaluated  as  a  weighted  sum  of  the 
quasi-steady  and  nonsteady  measurements,  and 
should  more  or  less  deviate  from  that  correspond¬ 
ing  to  the  classical  logarithmic  velocity  profile. 
For  the  study  of  the  nonsteady  regime,  frontal 
:heory  might  be  expected  to  play  an  essential 
role.  However,  such  a  study  will  not  be  attempted 
before  we  enter  the  second  phase  of  our  research 
program. 

FIELD  MEASUREMENTS 

Experimental  Considerations 

As  stated  earlier,  in  a  natural  field  over  a 
sufficiently  long  time,  the  entire  spectrum  of  eddy 
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scales  appears  and  hence,  theoretically,  the 
relevant  scale  encompasses  the  complete  spectrum. 
Such  longtime  averages  of  atmospheric  measure¬ 
ments  are  definitely  of  interest  in  the  study  of 
weekly  or  seasonal  variations  of  atmospheric 
phenomena  >md,  in  particular,  of  horizontal 
transport  processes.  On  the  other  hand,  over  a 
relatively  short  interval,  we  can  usually  identify 
two  distinct,  dominating  ranges  of  eddy  scales,  of 
which  the  smaller  one  alone  controls  the  vertical 
transport  process  in  the  earth’s  atmosphere. 
Thus,  inasmuch  as  we  restrict  our  study  to  the 
vertical  atmospheric  distribution,  the  long  period 
averages  can  be  ruled  out  as  improper.  The 
“relevant  scale”  for  the  vertical  distribution 
study  is  the  small  one;  the  corresponding  averaging 
time  is  the  short  periods;  i.e.,  of  the  size  of  meters 
and  period  of  10  seconds. 

Based  on  the  above  reasoning,  a  new  approach 
for  handling  this  problem  was  proposed:  separa¬ 
tion  of  the  transient  part  of  the  wind  from  the 
steady  part  and  a  study  of  their  characteristics 
as  separate  entities.  In  order  to  envision  more 
clearly  the  picture  involved,  one  may  consider 
the  atmospheric  boundary  layer  as  consisting  of 
three  distinct  cascading  layers,  as  follows: 

(1)  Laminar  sublayers,  at  the  surface  of  the  leaf 
and  stalk; 

(2)  Micro-meteorologic,  turbulent  boundary 
layer  (over  plant  top  and  canopy  flow  within 
plants) ; 

(3)  Synoptic,  turbulent  boundary  layer  of 
earth’s  atmosphere. 

These  three  layers  embrace  different  scales  of 
eddies  and  have  different  mean  velocity  gradients. 
The  roughness  of  leaf  and  stalk  surfaces  as  seer:  by 
the  laminar  sublayer  does  not  directly  influence 
the  mixing  mechanism  of  the  micro-turbulent 
layer.  On  the  other  hand,  the  fluctuating  velocity 
in  the  laminar  sublayer  is  directly  the  result  of  the 
eddies  of  the  micro-turbulent  layer.  In  the  same 
manner  the  roughness  of  plants  as  seen  by  the 
micro-turbulent  layer  plays  no  direct  role  in  the 
mixing  mechanism  of  the  larger  synoptic  atmos¬ 
pheric  boundary  layer.  However,  the  large  eddies 
of  the  synoptic  layer  are  responsible  for  the  gusti¬ 
ness  of  the  micro-turbulent  layer.  This  is  pre¬ 
cisely  why  the  micro-turbulent  layer  is  never 
steady.  The  turbulence  in  the  synoptic  layer  is 
controlled  essentially  by  the  roughness  scale  of 
much  larger  magnitude,  such  as  terrestrial  irregu¬ 
larities  (hills  ana  valleys),  other  disturbances  gen¬ 
erated  bv  temperature  gradients  (buoyant  mo¬ 
tion),  and  continental  pressure  fields 

The  small-scale  eddies  in  the  micro-turbulent 
layer  are  now  solely  responsible  for  the  vertical 
diffusion  of  carbon  dioxide,  vapor,  and  heat  around 
plants.  This  is  the  laver  that  is  to  be  studied  in 
detail,  in  order  to  understand  the  mechanism  of 
lower  atmospheric  diffusion.  As  mentioned  ear¬ 
lier,  this  layer  is  influenced  by  the  large  eddies  of 


the  outer  synoptic  boundary  layer,  and  therefore 
the  state  of  flow  varies  constantly  between  the 
transient  and  the  quasi-steady.  During  the  tran¬ 
sient  stage,  the  wind  shear  is  not  constant  over  the 
height,  so  that  the  logarithmic  wind  profile  does 
not  exist.  However,  with  eddies  of  sufficient  size, 
a  period  of  quasi-steady  state  exists  and  this  then 
enables  one  to  measure  the  steady-state  flow- 
parameters. 

As  a  first  step  toward  the  understanding  of  a 
complex  problem,  this  report  is  limited  to  the  study 
of  the  quasi-steady  state. 

The  wheatfield  selected  for  the  present  study  is 
situated  on  top  of  Mount  Pleasant,  5  miles  east  of 
Ithaca,  N.Y.  Measurements  of  wind  profiles 
were  taken  during  July  1960,  starting  when  the 
heads  of  wheat  were  level  with  the  top  leaves,  and 
continuing  up  to  the  time  when  the  heads  were 
30  cm.  above  the  top  leaves.  The  average  peak- 
plant-height  ranged  from  95  to  130  cm.  Average 
plant  density  was  60  stalks  per  square  foot. 

The  instrument  mast  was  located  near  the  center 
of  a  200-  by  200-meter  field.  Small-cup  anemome¬ 
ters,  9  cm.  in  diameter,  were  used.  They  were 
seated  at  10,  20,  40,  80,  and  160  cm.  above  the 
peak  level  of  crop.  Counters  for  all  anemometers 
were  read  simultaneously  at  every  30-second  in¬ 
terval.  Attempts  were  made  to  select  the  wind- 
profile  record  only  when  a  quasi-steady  wind  ex¬ 
isted.  That  is,  a  wind  record  was  selected  when 
it  was  preceded  and  followed  by  a  30-second  wind 
record  with  approximately  the  same  average  wind 
velocity  as  measured  by  the  top  anemometer. 

From  the  experimental  record  it  was  evident 
that  for  a  considerable  part  of  the  time  the  wind 
profiles  over  the  crop  were  transient. 

Five  representative  steady  wind-profiles  of  dif¬ 
ferent  mean-velocity  ranges  are  plotted  in  figure  6. 
One  important  characteristic  is  that  the  mean  ve¬ 
locity  gradient  at  the  peak  level  of  the  crop  changes 
only  slightly  with  increase  in  mean  velocity.  In 
addition,  the  characteristics  of  these  profiles  are 
insensitive  to  thermal  conditions  and  stages  of 
plant  growth. 

Logarithmic  Wind  Profiles 

Steady  wind  velocity  profiles  over  a  crop  can  be 
expressed  by 

U(Z,  (7o)=^P  log  Z>h  (32) 

where 

k—0A2,  von  K&rm&n’s  constant; 

Z= height  above  ground; 

h  =  peak  level  of  crop; 

f’o^inean  velocity  outside  of  boundary  layer 
(at  10-meter  height  based  on  logarith¬ 
mic  profile); 
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£/*=  In  (10)  •  friction  velocity; 


d— zero-plane  displacement; 
Zq= rough  ness  length; 

.  T=shear  stress; 
p=mass  density. 


U d,  and  Z0  are  to  be  determined  from  the  meas¬ 
ured  wind  profiles,  and  they  are  found  to  be  de¬ 
pendent  on  U0.  The  value  of  d  is  obtained  from 
equation  33,  where  the  mean  velocity  gradients  at 
two  different  heights  are  obtained  from  measure¬ 
ments. 


U0  is  computed,  after  d  is  known,  as  follows: 


t  Ui-U* 
(jW) 
(Z,-d) 


(34) 


l-h 


Z0  is  computed  after  d  and  L\  are  found. 

(35) 

When  the  three  parameters  have  been  computed 
for  each  experimental  wind  profile,  the  value  of  U0 
(at  10  meters  height)  can  be  computed  from 
equation  32. 

Values  of  Z o,  £/*,  and  ^  are  plotted  in  figure  7  as 

a  function  of  U0.  Both  ^  and  U„  are  found  to  be 

insensitive  to  free-stream  velocity.  On  the  other 
hand,  Z0  decreases  rapidly  with  increase  of  U0  and 
remains  constant  when  U0  increases  beyond  8 
meters/sec.  The  drag  coefficient 

e»-iinao)|.^-2(£)’ 

2 

is  plotted  as  a  function  of  U0  in  figure  7.  The 
reduction  of  CD  with  increasing  U0  is  evident. 

Canopy  Flow 

The  basic  differential  equation  for  canopy  flow 
(35)  is 

(36) 

where 


t(Z)  =eddy  viscosity,  a  function  of  Z; 
p= dynamic  viscosity,  a  constant; 

A(Z)  =  total  vertical  plane  area  of  plants/unit 
volume,  a  function  of  Z; 

Ti,=drag  coefficient  of  plant; 

p=mas8  density  of  air; 

Uh= mean  velocity  at  height,  h. 

Within  the  flow  range  considered,  C'D  can  be 
taken  as  independent  of  Reynolds  number. 
Equation  36  can  be  rearranged  as 


U Ji*A (2)C dpU*/2  ~ 

W  «(Z)+m 


d«(Z)  dU 
dZ  dZ 


^S(Z)  y  (37) 


where  77  and  Z  are  U  and  Z  normalized  by  U„ 
and  h,  respectively.  We  consider  S(Z)  to  be  the 
parameter  (shape  function)  of  the  differential 
equation.  For  wheat,  A(Z)  can  be  considered  as 
independent  of  (2),  because  the  area  density  of  the 
plant  is  practically  constant  with  height.  The 
eddy  viscosity  t(Z)  may  be  taken  as  being  a 
constant  for  the  top  quarter  of  Z,  for  the  reason 
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that  the  energy-carrying  eddies  are  not  expected 
to  be  modified  extensively  as  they  penetrate  into 
the  plants.  Hence  S(Z)  may  be  taken  as  a  con¬ 
stant  1§. 

With  boundary  conditions 


U\jm0—0 


equation  37  is  unique  and  can  be  solved  for  each 
value  of  S.  The  canopy  wind  profiles  are  plotted 
in  figure  8  for  various  shape  functions  from  10  to 


W 

;,000.  The  velocity  gradient  at  h,  — — 

dZ 


is 


plotted  in  figure  9. 

The  velocity  gradient  of  the  wind  profiles  at  h 
must  match  the  velocity  gradient  at  the  edge  of 
the  canopy  flow_(fig.  9).  From  this  condition  the 
shape  function,  S,  as  a  function  of  outer  boundary- 
layer  velocity,  U0,  can  be  found  and  is  plotted  in 
figure  7.  This  completes  the  tie  between  the 
wind  profile  and  canopy  flow'  once  the  steady  wind, 
U0,  is  given. 

The  shape  function,  5,  drops  rapidly  with  in¬ 
crease  of  and  tends  to  reach  a  terminal  value 
above  a  7-meter/sec.  windspeed.  In  figure  S,  one 
notes_that  the  canopy  profiles  are  more  curved  for 
high  S  and  low  U„,  and  approach  asymptotically  a 


Relative  Height  (Z) 


Fnu  re  8. — Canopy  wind  profiles  for  various  constant 
shape  functions.  .S. 


Jfii 

tilh 


Figure  9. —  Velocity  gradients  at  peak -plant-height  as 
a  function  of  shape  function,  S. 

straight  line  for  low  S  and  large  U0.  S  varies  di¬ 
rectly  as  CV ,  and  inversely  as  «(Z);  the  reduction 
of  S  at  large  U0  means  a  reduction  in  fD'  and  an 
increase  in  t(Z).  This  in  turn  means  that  the 
flow  is  more  turbulent.  We  see  therefore  that 
the  field  measurements  exhibit  trends  similar  to 
those  anticipated  from  the  theory. 

Accurate  measurement  of  canopy  flow  has  not 
yet  been  made,  but  from  some  preliminary  meas¬ 
urements  there  seems  to  be  a  tendency  for  the 
canopy  wind  profile  to  cross  from  a  lower  S’  value 
to  a  higher  S’  value  curve,  as  Z  goes  below  O.S. 
This  would  mean  that  e(Z)  may  not  be  a  constant, 
but  should  be  reduced  as  Z  is  lowered.  Thus  far, 
experimental  data  are  not  sufficient  to  evaluate 

<(Z). 

From  the  known  value  of  the  S  factor,  for  a 
given  f.\  and  plant  density,  the  value  of  eddy 
viscosity  can  be  computed.  For  the  range  of  data 
presented,  the  ratio  of  eddy  viscosity  to  dynamic 
viscosity  ranges  from  1  X  103  to  3X  104.  The  eddy 
viscosity  at  peak  plant  height,  <»,  can  be  com¬ 
puted  from  the  logarithmic  equation,  equation  32, 

tH=Plc(h-d)CJ\n  (10).  (38) 

The  ratio  of  eddy  viscosity  to  dynamic  viscosity 
obtained  by  this  method  is  considered  to  be  more 
accurate,  and  it  was  found  to  range  from  0.9  X  HV 
to  2.2X103.  The  values  at  high  range  do  not 
agree  well  with  those  computed  by  canopy  flow. 
The  error  is  probably  due  to  the  assumption  of 
constant  t(Z).  More  experimental  measurements 
and  analytical  study  on  <(Z)  will  be  conducted. 


A  STUDY  OF  ATMOSPHERIC  TURBULENCE  AND  CANOPY  FLOW 


By  H.  S.  Tan  and  S.  C.  Lino 


In  the  present  study  the  following  two  aspects 
are  stressed:  (1)  To  furnish  the  “relevant  scale” 
concept  with  a  more  rigorous  foundation  and  (2) 
to  improve  the  canopy  flow  model  bv  postulating 
the  linear  growth  of  diffusivity  rather  than  the 
logarithmic  velocity  profile  as  the  basic  property 
of  turbulent  wall  flow.  It  is  gratifying  to  find 
that  the  predicted  velocity  profile  curvature  re¬ 
versal  deep  in  the  canopy  flow  was  indeed  ob¬ 
served  and  pointed  out  repeatedly  by  several 
investigators  (22). 

It  might  not  be  out  of  place  to  mention  at  this 
point  that  until  now  we  have  not  approached  the 
diffusion  problem  with  source  distribution,  the 
photosynthesis  problem,  and  the  thermo-insta- 
nility  problem  as  seriously  as  we  would  like  to. 
The  main  reason  behind  this  is  the  heavy'  de¬ 
pendence  of  these  phenomena  on  the  atmospheric 
turbulence,  of  which  little  is  definitely  known  and 
much  has  yet  to  be  clarified  before  we  can  satis¬ 
factorily  formulate  a  working  model  of  atmospheric 
transport. 

An  introductory  discussion  of  turbulent  dif¬ 
fusivity  in  shear  flow  is  provided  by  Dr.  Spence. 

Our  study  so  far  has  been  restricted  to  the 
quasi-steady  part  of  the  mean  motion  (time  mean). 
A  tentative  shear  wave  model  for  treating  the 
nonsteady  problem  has  been  suggested  by  Dr. 
Ling,  and  included. 

TURBULENT  TRANSPORT  IN  LOWER 
ATMOSPHERE 

Equation  of  Turbulent  Transport 

The  transport  of  any  physical  quantity  in  the 
atmosphere  is  governed  by  the  following  basic 
differential  equation: 

ft  =  £t+qW=V.(KV8)+S(e)  (39) 

in  which  B  denotes  the  transported  quantity,  q 
the  wind  velocity,  K  the  coefficient  of  diffusivity, 
and  S  the  source  function.  Physically  the  equa¬ 
tion  states  that: 

(a)  The  total  derivative  of  quantity  B  is 
given  by  the  sum  of  divergence  of  flux 
and  rate  of  production  of  B; 

(b)  The  flux  of  B  is  in  turn  proportional  to  the 
gradient  of  B\ 

673752  a  -63 - 3 


(c)  The  flux  per  unit  gradient  is  named  the 
coefficient  of  diffusivity. 

In  connection  with  this  transport  equation,  the 
following  two  studies  are  of  particular  interest: 

(a)  Simplification  of  the  equation  to  corre¬ 
spond  with  specific  physical  problems; 

(b)  Effect  of  turbulence  on  diffusivity,  K, 
and  hence  on  its  solution  B. 

For  lower  atmospheric  study  the  following 
simplifications  are  usually  obtained. 

1.  Vertical  Gratification  and  horizontal  wind 

V=k  q—tfj:  so  q  ■  V0* 0 

and  equation  39  becomes 

(40> 

which,  upon  putting  *9=0,  forms  the  basic  equa¬ 
tion  for  usual  nonsteady,  micro-meteorology 
study: 

2.  Steadiness  condition, 

ot 

Equation  39  further  reduces  to 

Hk£)+S(I)=°  <41) 

Equation  41  actually  forms  the  basis  for  the 
first  two  sections  of  this  report  (pp.  3  and  7) 
and  for  the  present  investigation.  Upon  putting 
S= 0,  equation  41  simplifies  to 

al«)=0;  (4U) 

i.e.,  the  well-known  “constant  flux”  equation  for 
steady  earth  surface  layer. 

For  the  present  report,  B  will  denote  either 
carbon  dioxide  concentration  in  photosynthesis 
study  or  horizontal  windspeed,  u,  in  velocity  dis¬ 
tribution  study. 

It  is  clear  from  the  above  that  our  study 
naturally  separates  into  two  regions  by  the  crop 
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height,  z=h.  For  the  outer  flow  (2>A),  o-O, 
the  governing  equations  for  both  (COjj  and  u  are 
given  by  equations  40a  and  41a;  i.e.,  the  well- 
known  atmo^oheric  diffusion  equation.  For 
quasi-steady  cas# ,  we  have  simply: 


For  the  inner  flow  (2<A),  on  the  other  hand, 
SV  0.  In  this  case  it  has  been  shown  that  for 
quasi-steady  flow  we  have  (35): 

[CO,]  study:  ^  (if  ™)+(Xz)9+g(z)=0;  (43) 

u  study:  —  (if  ^^+S(2)tf,=0.  (44) 

Turbulence  and  Diffumity 

The  role  played  by  turbulence  in  atmospheric 
transport  takes  on  the  form  of  coefficient  of 
diffusivity,  K.  Prandtl  first  introduced  his  fa¬ 
mous  mixing  length  theory  for  the  mechanism  of 
turbulence  transport  and  showed  that  if  is 
actually  proportional  to  his  mixing  length,  /.  By 
further  postulating  that  in  the  neighborhood  of 
a  solid  wall,  mixing  length,  l  (and,  consequently, 
diffusivity,  K)  grows  linearly  with  distance  z  from 
the  wall,  he  was  able  to  deduce  from  equation  42 
that  the  velocity  profile  for  turbulent  wall  flow 
must  follow  a  logarithmic  law;  i.e., 

U=T (f.)  <«> 

This  velocity  profile  has  indeed  been  confirmed 
by  many  observations. 

In  the  last  20  years,  however,  there  has  been 
a  noted  trend  among  investigators  toward  taking 
a  diametrically  opposite  view  in  discarding  the 
mixing-length  hypothesis  completely  and  accept¬ 
ing  the  logarithmic  law  as  the  basis  for  theoretical 
development,  from  which  the  “austausch  coeffi¬ 
cient”  K  can  be  shown  to  grow  linearly  with  z. 
So  far  as  the  outer  flow  is  concerned,  it  makes 
no  difference  which  point  of  view  to  take.  •  How¬ 
ever,  as  soon  as  we  enter  into  the  inner  flow 
region,  where  there  are  distributed  sources,  both 
approaches  face  impossible  difficulty.  On  the 
other  hand,  by  extending  Prandtl’s  linear  growth 
hypothesis  instead  to  the  canopy  flow  diffusivity, 
i.e.,  /f~2,  equation  44  becomes: 

<44»> 

So  we  were  able  to  compute  the  theoretical  velocity 
profile  and  predict  a  curvature  reversal  deep  in 


the  cro*\  which  has  been  amply  verified  by  various 
field  measurements. 

This  fact  is  considered  to  be  of  great  significance 
toward  justifying  Prandtl’s  mixing  length  concept 
in  general,1  and  pointing  to  the  basic  correctness 
of  our  theoretical  canopy  flow  model  in  particular, 
i.e.,  turbulent  flow  with  linearly  growing  diffusiv¬ 
ity,  If. 

Theory  of  Average  and  Its  Applications 
1.  Theorem* 

The  following  definitions  and  theorems  are  in¬ 
tended  for  clarifying  some  basic  questions  con¬ 
cerning  the  process  of  ensemble,  space,  and  time 
averages  employed  in  the  atmospheric  turbulence 
study. 

Definition. — The  mean  value  of  a  bounded  fluctu¬ 
ating  function  /  is  defined  by: 


*  1  /*<+*• 
S(»l'T=5Tjt_r  MW- 

Theorem  1. — It  is  evident  that 

(a)  >«)lr-H>  =./«); 

(b)  ,/(<)|  r-*.=const.; 

(c)  /(f)|r  =/n(f,r). 

(47) 

Theorem  2. — The  mean  of  a  simple  periodic 
function  retains  its  periodicity,  T0,  but  the  ampli¬ 
tude  falls  off  according  to  (T/To)-1.  It  may  be 
noted  that  there  exists  a  set  of  eigen-functions: 

/(O|r-»r0=O- 

(48) 

Definition. — The  spectrum  of  a  fluctuating  func¬ 
tion/  is  defined  by: 

JM\T=joT  .me'-dt^fi^D) 

(49) 

/(*)U=J]  j\x)etudx=fik,X). 

(49a) 

Theorem  3. — It  is  evident  that  /(«,<*> )=/(w), 
/(jkco)=/(it);  i.e.,  the  Fourier  Spectra. 

Definition. — We  define  the  two  following  time 
and  space  correlation  functions: 

ft(r)==L  u(x,  t)u(x,t+r);  (50) 

u '* 


*  Indeed  it  would  be  interesting  to  introduce  a  parametric 
representation  of  the  “Austausch"  coefficient  as  follows: 
A’~r*  where  a  is  a  function  of  S.  a—* >1  as  S— d).  a(S)  can 
be  determined  by  correlating  the  canopy  flow  velocity 
profile  with  varying  crop  density. 
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R(i)=^jHx,t)u(x+i,t);  (51) 


where  (“)  denotes  ensemble  average. 

Theorem  4- — Ergodic  Principle:  If  a  turbulent 
field  has  a  steady  time  mean  and  a  uniform  space 
mean,  then  the  field  is  ergodic ;  that  is,  the  ensemble 
average  can  be  replaced  by  either  time  or  space 
average. 

In  the  present  case  this  gives 


/?(r)=Lim=Li  f  u(x,t)u(x,  t+r)dt ;  (50a) 

T-*m  U '*  1  JO 

R(()=Um  JL>i  f  u(x,t)u(x+(,t)dt.  (51a) 

r->.  u  * 1  Jo 

Infinite  extent  of  T  or  X  cannot  be  reached, 
and  so  a  real-steady,  uniform  mean  does  not 
exist.  This  state  can,  however,  be  approximated 
by  a  quasi-steady  and  quasi-uniform  mean  over 
T  and  7  denotes  the  rr  \imum  averaging 
period  or  length  of  record  and  X  the  actual  field 
extent. 

For  a  turbulent  field  with  quasi-steady,  quasi- 
uniform  mean,  quasi-ergodic  hypothesis  permits 
the  following  replacement: 

yn-^joT.f^.t)dt  (52) 


Lim  RT  ((D=B(t,  {).  This  means  Rt  t(T)  is 
r->« 

a  faster  decreasing  function  than  R,  (  with  increas¬ 
ing  r,  {. 

Hence,  by  the  Birkhoff  Ergodic  Theorem  {12), 
the  applicability  of  a  quasi-ergodic  hypothesis  is 
assured. 

Theorem  7. — Eddy  scale  L  is  defined  by  (27): 

L= j‘  R((T)d(=L({T)  (56) 

=arj*  Rf{T)dT—alj,{T)  (56a) 

from  theorem  6.  Thus,  the  recognizable  eddy 
scale  grows  with  the  interval  of  average  T. 
Theorem  8. — In  equation  56a,  we  can  put 

a  =  U.  (57) 

By  field  measurements,  Cramer  (4)  established  the 
following  remarkable  relationship, 


Rt(T)=R((T),1=Ut,  (58) 

based  on  a  20-minute  sampling  interval  ( T)  with 
U  the  20-minute  mean  windspeed.  These  cor¬ 
relation  curves  fit  approximately  the  Kolmo- 
goroff  Law  for  inertial  subrange;  i.e., 


where  T  is  the  finite  averaging  period,  xH—x(tn), 
and  f,_i>7’.  This  gives — 

Theorem  5. — 

5(r)  =  =L  TTrn  2  f  u(xH,  t)u(XH  t+ r)dt  =  Rr(T)  J 
U  *  •'»  ■*  B-l  Jo 

(53) 

^({)=4iA7fS  f  u{xn,t)u(xn+l,t)dt=R({T). 

U  *  1  B-l  JO 

(54) 

Theorem  6. — For  a  bounded,  fluctuating  func¬ 
tion  f(x,  t)  with  quasi-steady,  quasi-homogeneous 
mean,  we  postulate  the  following  similarity  law: 

RriD^RdDl  (55) 

i.e.,  as  functions  of  averaging  period,  T,  the  space 
and  time  correlation  functions  behave  in  similar 
manner. 

The  general  behavior  of  R,  and  /fj,  as  functions 
of  T  in  the  lower  atmospheric  turbulence,  is 
shown  in  reference  (24:  hg.  i.2.1,  p.  4<J)  and 
reference  (3:  Jig.  8.4-15 ,  p.  201),  respectively. 
From  the  figures  given  in  the  references,  it  is 
evident  that  both  Rr  and  /?{  are  monotonically 
increasing  functions  of  T,  and,  indeed, 


(59) 

and  the  agreement  improves  with  increasing 
turbulence  level,  indicating  that  isotropy  extends 
further  into  lower  frequency  range  with  higher 
turbulence  intensity. 

2.  Application  to  atmocpherie  turbulence 
analytic 

An  atmospheric  wind  record  exhibits  random 
fluctuations  containing  many  freauencies,  and 
usually  a  few  dominant  ones.  The  maximum 
averaging  period,  T,  is  limited  from  above  by  the 
length  of  the  record  available. 

An  instrument  has  finite  resolution  and,  hence, 
can  only  read  a  corresponding  average  with 
7>0.  The  minimum  averaging  period,  T,  is  thus 
limited  from  below  by  the  instrument  sensitivity. 

An  ideal  spectrum  \fu^(u>)  corresponds  to 
T-»®,  T-d).  Now  what  are  the  effects  of 
finiteness  of  7  and  T  on  a  spectrum? 

Theorem  5  shows  that  the  maximum  recog¬ 
nizable  eddy  scale  is  limited  by  T.  [Indeed 
L(T)  is  a  monotonically  increasing  function  of  T.\ 
Thus,  7  effectively  offer  a  low  frequency  (low- 
wave  number)  cutoff,  ’  the  sense  that  Deyond 
this  point  the  spectrum  function  defined  by 
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equation  39  only  holds  formally  and  loses  its 
physical  significance. 

'theorem  2  and  equation  49  show-  that  the 
multiplying  factor  for  averaging  over  period  T 
of  a  periodic  function  with  periodicity  T0  varies 
as  (i/To)-1.  Thus,  because  of  the  inevitable 
averaging  character  of  an  instrument,  the  spec¬ 
trum  of  higher  frequency  (7,0<7T)  or  small- 
scale  fluctuating  components  are  effectively  sup¬ 
pressed.  In  other  words,  this  averaging  process 
removes  a  large  part  of  the  time  dependence 
due  to  small-scale  eddies,  but  retains  the  time 
dependence  of  large-stale  eddies.  The  T  of 
instrument  sensitivity,  therefore,  imposes  the  high 
frequency  (high-wave  number)  cutoff. 

In  relation  to  the  relevant  eddy  scale,  L(T), 
for  vertical  atmospheric  turbulent  transport  proc¬ 
esses,  the  smaller  scale  eddies  play  the  role  of  turbu¬ 
lent  mixing  (austausch),  while  larger  eddies  appear 
as  nonsteady  wind. 

If  there  were  no  dominant  scales  of  atmospheric 
eddies,  then  the  word  “mean”  would  have  little 
meaning.  In  the  atmospheric  surface  layer, 
however,  there  are  usually  two  natural  dominant 
scales  of  eddies:  the  small-scale  ones  produced  by 
the  roughness  elements,  such  as  plant  leaf  and 
stalk;  and  the  large-scale  ones  produced  by  the 
terrain  irregularities,  like  hill  and  dale.  The 
existence  of  these  definitely  separated  ranges  of 
eddy  scale  has  been  demonstrated  by  field  record ; 
it  gives  a  unique  meaning  to  the  word  “mean,” 
ana,  moreover,  it  introduces  in  a  natural  way  the 
concept  of  relevant  scale  and  proper  averaging 
period. 

ANALYSIS  OF  ATMOSPHERIC  WIND 
IN  THE  VICINITY  OF  PLANTS 

Background  Information  of  This  Study 

The  atmospheric  wind  consists  of  a  cascade  of 
turbulent  eddies,  varying  in  size  from  a  few  meters 
at  the  plant  to  macroscale  at  the  fringe  of  the 
atmosphere  (36).  The  wind  in  the  vicinity  of  a 
crop  is  influenced  by  all  scales  of  atmospheric 
disturbances.  Thus,  strictly  speaking,  no  station¬ 
ary  time  average  exists.  However,  if  we  limit 
our  interest  fo  i  finite  scale  and  period,  a  quasi- 
steady  average  can  usually  be  found. 

The  basic  turbulent  mechanism  in  the  vicinity 
of  a  crop  is  controlled  primarily  by  the  roughness 
effect  of  the  crop  alone.  In  addition,  this  mech¬ 
anism  is  maintained  by  u  driving  wind,  l ’0,  that, 
within  a  given  period  of  interest,  can  be  repre¬ 
sented  by  a  mean  wind  and  Fourier  spectrum  of 
oscillating  winds.  This  driving  wind  is  a  function 
of  all  external  influences  other  than  that  of  the 
crop.  With  this  in  view,  it  is  possible  to  treat 
the  turbulent,  atmospheric  boundary  layer  as 
three  intermixed  states  of  flow,  namely:  (l! 
steady;  (2)  quasi-steady  ;  and  (3)  nonsteady. 


The  mean  driving  wind  contributes  to  the  local 
steady  state,  while  the  short-wave  spectrum  of  the 
fluctuating  wind  constitutes  its  turbulent  mech¬ 
anism.  A  sustained  local  wind  or  a  long  wave  of 
1,000  meters  can  be  taken  as  quasi-steady  flow. 
A  wavelength  in  the  range  of  40  to  1,000  me^e*^  is 
considered  nonsteady.  The  oscillating  compo¬ 
nents  in  this  range  propagate  downward.  Very 
little  is  known  about  nonsteady  turbulent-bound¬ 
ary -ln>er  flow.  Indeed,  the  types  ot  nonsteady 
Ijoundarv  layers  investigated  are  either  laminar- 
boundary  layers  ( 8 ,  17,  27,  38)  or  turbulent¬ 
boundary  layers  (11)  driven  by  a  nonturbulent, 
oscillatory  flow  field  (potential  flow).  A  bound¬ 
ary  layer  driven  by  a  tut  Indent  shear  flow  remains 
a  subject  yet  to  be  studied. 

The  most  important  feature  that  governs  the 
turbulent  transport  in  the  vicinity  of  a  crop  is  the 
steadv-state-velocity  profile.  To  obtain  these  pro¬ 
files  by  field  measurement,  the  sampling  period 
required  should  be  greater  than  the  relaxation 
period  of  the  boundary  layer.  The  relaxation 
period  will  be  defined  in  the  section  on  shear-wave 
theory.  For  flow  over  a  crop,  this  relaxation  time 
is  approximately  10  seconds.  An  averaging  period 
of  30  seconds  should  smooth  out  the  high-fre¬ 
quency  waves  due  to  turbulent  eddies.  If  the 
averaging  period  is  extended  to  cover  the  maxi¬ 
mum  range  of  the  nonsteady  wind,  the  time  de¬ 
pendent  factor  of  this  component  can  again  be 
eliminated.  However,  it  can  be  shown  from  field 
measurements  (fig.  10)  that  the  intensity  of  non¬ 
steady  as  well  as  quasi-steady  winds  is,  in  general, 
beyond  the  linear  perturbation  range.  Conse¬ 
quently,  the  result  of  a  long  averaging  period  is  to 
include  nonlinear  effects  in  the  steady-state 
average. 

So  far,  the  short-period  averaging  seems  to  be 
most  logical  for  this  study,  but  this  must  be 
accompanied  by  a  selective  sampling  technique 
to  eliminate  the  nonsteady  components  as  we’,1 
as  the  nonlinear  effects  due  to  extreme  variation 
of  wind  velocities.  The  method  used  in  this 
report  embraces  the  method  of  selecting  samples 
from  periods  of  30-second  duration  preceded  and 
followed  by  approximately  the  same  time  average. 
Field  data,  obtained  bv  this  means,  were  found  to 
be  quite  consistent.  This,  in  a  way,  verified  the 
basic  concept  of  selective  short-time  average. 

Wind-profile  measurements  over  a  wheat  field 
und  a  cornfield  were  conducted  with  multiple 
small-cup  anemometers  of  9 -cm.  diameter,  in 
conjunction  with  Hasting  thermocouple-anemom¬ 
eters.3  The  cup  anemometers  were  appropri¬ 
ately  mounted  at  various  heights  within  and 
above  the  crop.  Thermocouple  anemometers 
were  used  within  the  crop  where  the  wiudspeed 
was  considered  too  low  for  cup  anemometers. 

•  The  use  of  this  or  other  patented  equipment  in  this 
study  does  not  imply  approval  of  the  product  to  the 
exclusion  of  others  tnat  may  also  be  suitable. 


»  Figure  10. — Simultaneous  horizontal-wind  recording  at  two  different  elevations  above  a  corn  crop. 

Simultaneous,  integrated  readings  from  these  of  Agriculture,  under  the  direction  of  E.  R.  Lemon, 
,  anemometers  were  recorded  at  30-second  intervals.  ARS,  USDA.  A  cons  lerable  amount  of  field 

,  The  average  wind  profiles  were  obtained  from  data  had  been  compiled  during  the  summer 

|  these  measurements.  The  selected  samples  rep-  months  of  1960,  but  only  a  small  part  of  the 

resenting  quasi-steady  wind  profiles  for  a  wheat  information  was  analyzed  and  presented  herein,  as 
,  field  and  a  cornfield  are  shown  in  figures  1 1  and  12,  a  consequence  of  limited  time  and  effort. 

>  respectively.  The  range  of  wind  profi,  pre-  The  general  theory  and  experimental  results 

sented  here  is  limited  to  the  most  common  pre-  of  steady  state  and  quasi-steady  state  winds  in  the 
j  vailing  winds  encountered  in  a  typical  field.  vicinity  of  a  crop  are  presented  in  the  next 

,  Data  for  extremely  low  and  high  winds  are  in-  subsection. 

k  complete,  and,  therefore,  not  presented. 

The  wheatfield  was  situated  on  top  of  Mount  Boundary-Layer  Profiles 

Pleasant  unobstructed  on  all  sides,  whereas  the  , 

I  cornfield  was  situated  in  Ellis  Hollow.  (See  L  Logarithmic  wind- profile  over  a  crop 

,  last  section,  p.  37.)  The  wind  profile  character-  The  logarithmic-boundary-layer  theory  of 

|  istics,  by  virtue  of  the  sampling  technique  em-  Prandtl  {23)  applies  very  well  to  the  steady- 

ployed,  are  expected  to  be  independent  of  the  boundary-layer  now  above  a  crop.  In  the  original 
geographical  terrain.  theoretical  model,  it  was  assumed  that  the  mixing 

|  Most  of  the  field  data  and  the  instrument ation  length  or  the  equivalent  eddy  viscosity  of  turbulent 

work  used  in  this  report  were  carried  out  with  the  mixing  grows  linearly  with  the  height.  This 

assistance  of  J.  Stoller,  New  York  State  College  assumption  is  now  extended  all  the  way  from 
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Figcre  11. — Wind  profiles  over  a  wheat  crop,  I960. 
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Figcre  12. — Wind  profiles  over  a  corn  crop.  1960. 

and  incorporating  equations  61  and  62  in  equation 
60,  we  obtain  the  well-known  logarithniic- 
boundary-layer  equation  through  integration  of 
equation  60. 


within  the  crop  through  the  inner  boundary  layer 
up  into  the  outer  turbulent-flow  layer  (21).  By 
virtue  of  our  stratification  hypothesis;  i.e.,  van¬ 
ishing  pressure  gradient^  >  we  have 

-=0- 
dZ  ’ 

i.e.,  the  flow  is  maintained  through  the  shearing 
action  of  eddv  viscosity.  Now  this  leads  to 


r/,ar  t  . 

r  =  t(Z;  -  7,=constant. 
oZ 


(60) 


L(Z,U0)  1.  Z  d(Li a)  j 

T .(Gi  '^'W  z>h  m 


where 


&  =  0.42  Von  Karman’s  constant; 

Z=  height  above  ground  level; 
f mean  wind  velocity  at  10-meter  reference 
level ; 

h  —  peak  plant  height,  defined  as  1.1  times 
the  mean  plant  height,  A; 

<i  —  zero-plane  displacement,  a  function  off’,; 
Z„= roughness  length,  a  function  of  f’„; 
p—  mass  density  of  air. 


Introducing  the  friction  velocity  as 

r*s\  p’ 

the  eddy  viscosity  «(Z),  as 

t(Z)  =  pl\k(Z-d),  (62) 


The  constants,  f  *,  Z„,  and  d  are  called  the  char¬ 
acteristic  parameters  of  a  steadv-state-wind  profile. 

In  order  to  provide  a  convenient  reference  level 
in  place  of  a  boundary  thickness,  the  conventional 
height  of  10  meters,  Zl0,  was  chosen  as  the  nominal 
limit  of  the  surface  layer.  Ratios  formed  with 
Z ,0  are  to  be  taken  as  of  approximate  order  of 
magnitudes. 
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The  three  wind-profile  parameters,  d,  ,  and  Zo, 
can  be  determined  from  the  measured  wind  pro¬ 
files,  as  shown  in  figures  11  and  12,  and  they  are 
obtained  through  the  following  equations: 


and 


ran  7 
LazJ, 

i  i 

oi  O' 

Z, 

1 

fatn 

Lazi 

ratr 

Laz. 

1 

n  - 

i\=i  Ll 

-u2 

(Z.-rf) 


in 


7  d 

z'»~  e*Ui/v*  ' 


164) 


(65) 


(66) 


These  equations  are  derived  directly  from  the 
basic  logarithmic  equation  63.  The  values  of  these 
parameters  as  a  function  of  reference  wind,  U0,  for 
a  wheatfield  and  a  cornfield  are  plotted  in  figures 
13  and  14,  respectively.  It  is  to  be  noted  that 
values  given  in  the  above  figures  are  for  a  typical 
wheatfield  and  a  cornfield  having  an  approximate 

fdant  density  of  60  plants  and  0.5  plant  per  square 
oot  of  ground,  respectively. 

The  value  of  kinematic-eddy  viscosity,  Km, 
known  also  as  the  coefficient  of  eddy  diffusivity, 
is  of  basic  interest  in  the  study  of  diffusion  of 
COj,  water  vapor,  and  heat.  Direct  measure¬ 
ment  of  this  coefficient  was  planned,  but  only 
crude  measurements  were  possible  at  the  time 
of  this  report.  In  general,  they  check  well  with 
those  obtained  indirectly  from  the  wind-profile 
parameters.  The  equation  for  eddy  diffusivity 
is,  from  equation  62, 


Km=^=U,(U0)k[Z-dm).  (67) 

p 


Values  of  t/*  and  d  can  be  obtained  from  figures 
13  and  14  for  a  given  U0  and  field  condition. 

The  drag  coefficient  defined  as 

(6S) 

2 

is  also  plotted  in  figures  13  and  14. 

Comparing  the  behavior  of  flow  parameters 
between  a  wheat  and  a  cornfield,  it  is  noticed 
that  the  zero-plane  displacement,  d,  decreases 
with  an  increase  of  U„;  this  characteristic  is 
controlled  primarily  by  the  mechanism  of  the 
canopy  flow.  The  values  of  friction  velocity, 
t  *,  are  of  about  the  same  magnitude  and  char¬ 
acteristic  for  both  the  wheatfield  and  the  cornfield 
This  indicates  that  the  characteristics  of  turbu¬ 
lence  above  either  crop  are  of  approximately  the 
same  scale  and  intensity.  The  value  of  roughness 


parameter,  Z>,  behaves  differently  for  the  two 
crops.  For  the  wheatfield,  Zo  tends  to  decrease 
with  U0,  whereas  the  reverse  is  true  for  the 
cornfield.  It  is  postulated  that  the  behavior  of 
Z0  depends  on  the  flexibility  of  the  plant  as  well 
as  plant  density.  The  wheat  plant  is  more 
flexible  and  tends  to  bend  and  “streamline” 
itself  under  a  high  witui,  but  the  -^ru  is  relatively 
stiff  and  does  not  “streamline”  itself  as  readily. 
The  difference  in  effective  plant  density,  or  leaf 
area  density  between  wheat  and  com,  exercises 
a  strong  influence  on  the  canopy  flow  that  in¬ 
directly  affects  the  characteristic  of  Zo.  The 
different  dependence  of  CD  on  U0  between  a 
wheatfield  and  a  cornfield  is  similar  to  the  varia¬ 
tion  of  Zo  with  U<>.  Therefore,  the  same  tentative 
explanation  regarding  plant  flexibility  is  applicable. 

2.  Canopy  (low  within  a  crop 

The  basic  analytical  model  for  a  canopy  flow 
has  been  mentioned  in  the  introductory  section 
of  this  report.  Basically,  it  is  assumed  that  the 
flow  within  the  crop  is  maintained  by  the  action 
of  eddy  viscosity,  and  that  the  mixing  length 
theory,  which  is  applied  to  the  outer  flow,  is 
also  valid  in  this  zone.  The  differential  equation 
governing  the  flow  is, 

^C'oA(Z)^  (69) 

where 

A{7)  =  total  vertical  plane  area  of  plants/unit 
volume,  a  function  of  Z;  and 

( ©=drag  coefficient  of  plant. 

The  term  on  the  right  represents  the  drag  force 
offered  by  leaves  ami  stalks  of  the  crop.  Within 
the  flow'  range  considered,  C'D  can  be  taken  as 
independent  of  Reynolds  number,  and  the  drag 
force  is  proportional  to  the  square  of  the  velocity. 

Introducing  equation  60  into  equation  69, 
we  have 

al["z)i]-r^<Z)ST-  (70> 

Let 


normalized  height  variable,  and 

normalized  wind,  whore  f  /(  —  wind  at 
peak-plant  height, 


and  introducing  these  new  variables  in  equation 
70;  we  have 


df'  ,2r  Uy)r'  P U 

=+*(Z)  -—^/dLvKZK  D  9-’ 
*  dZ  z 


dZ 


Z<1. 


(71) 
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Fh.i  re  13. — Parameters  for  steady-state  wind  profiles  for  winter  wheatfield. 
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In  order  to  define  statistically  the  value  of 
peak-plant  height,  some  sampling  data  on  the 

?robabUity  density  distribution  of  com  height, 
'(h),  were  taken  and  plotted  in  figure  15.  The 
peak-plant  height  has  a  value  of  about  1 .2  times 
the  mean  height  of  plants.  However,  instead 
of  taking  this  figure,  a  value  of  1.1  times  the  mean 
height  of  plants  is  taken  as  the  definition  of  h. 
It  was  assumed  that  this  definition  is  also  appli¬ 
cable  to  wheat. 

The  area  distribution  function,  A(Z),  can  be 
obtained  statistically  for  a  given  crop.  However, 
for  simplicity,  it  was  assumed  a  priori  that  the 
average  vertical  projection  area  of  a  single  plant 
(with  the  vertical  plane  making  an  angle  of  45° 
between  the  planes  of  maximum  and  minimum 
projected  area)  has  an  area  distribution  that  is 
uniform  from  the  root  of  a  plant  to  85  percent  of 
its  height  ;  from  that  point  on,  the  area  reduces 
uniformly  with  height  to  zero  at  the  tip.  By 
combining  this  standard  area  function  with  the 
plant  height  distribution  function,  the  statistical 


area  distribution  function,  A(Z),  is  obtained  and 
plotted  as  a  ratio  of  area  density,  A 0,  at  the 
midplant  height.  This  function  is  shown  with 
“x”  marks  in  figure  15.  The  area  distribution 
function  again,  for  the  sake  of  simplicity,  was 
approximated  bv  a  straight  line  shown  in  the 
above  figure.  This  simplified  area  function  gives 
zero  area  at  1.1  times  the  mean-plant  height, 
and  hence  this  height  was  used  for  the  definition 
of  peak-plant  height. 

If  we  introduce  the  linear  growth  of  eddy  vis¬ 
cosity  with  height, 

e(Z) =Zf/„ 

where  tk is  the  eddy  viscosity  at  height,  k,  and  with 
the  above  simplified  function  of  A(Z),  equation  71 
is  written  for  two  zones : 

for  0<Z<0.73 
A(Z)=A0,  constant 
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Figuu  15. — Plant  area  and  height  distribution  functions 
for  corn. 
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(73) 


where  5"  is  the  dimensionless  parameter  of  the 
differential  equation.  It  is  to  be  called  the  shape 
factor  of  canopy  flow,  namely, 


_  ph'U  iAqCd 
tk 


(74) 


Equations  72  and  73  are  nonlinear,  differential 
equations  of  the  second  order.  They  are,  how¬ 
ever,  unique  and  can  be  solved  by  numerical 
integration  for  the  following  boundary  conditions  : 

f/2-,  =  l; 

Uz.  o=0. 

Relaxation  method  (2S)  wras  used  to  obtain  the 
characteristic  canopv-wind  profiles  for  various 
values  of  shape  factors  The  _result,  plotted  in 
figure  16,  covers  the  range  of  S  from  5  to  1,000 
In  all  these  cases,  the  ground  shear  was 
neglected  in  comparison  to  the  leaf  drag.  For  S 
<5,  ground  effect  comes  into  play.  When  S-* 0, 
ground  shear  alone  controls  the  flow,  and  we 
obtain  the  normal  logarithmic  wind  profile. 
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F lot  re  1G. — Canopy  wind  profiles  for  various  shape 
factors  S. 


3.  Combination  of  canopy  and  logarithmic- 
wind  profile 

The  canopy  flow  and  the  logarithmic-profile  flow 
above  are  actually  one  continuous  system.  The 
matching  point  between  the  two  flow  fields  is  at 
the  peak-plant  height,  h,  where  the  wind  velocity 
and  the  velocity  gradient  must  match.  This 
means  that,  for  a  given  logarithmic-wind  profile, 
there  will  be  a  canopy  profile  of  a  given  shape 
factor  that  will  match  it  at  h.  The  value  of  «»  in 
equation  74  must  also  match  the  value  of  «(Z)» 
determined  by  equation  62  of  the  logarithmic-wind 
profile 

dC 

The  velocity  gradient  — _  for  the  canopy 

dZ  z-  1 

flow  is  plotted  in  figure  17  as  a  function  of  shape 
factor  .S'.  From  that  plot,  the  matching  shape 
factors  tor  the  wind  profiles,  shown  in  figures  1 1 
and  12,  are  plotted  in  figures  13  and  14,  respec¬ 
tively.  The  trend  of  shape  factor  for  the  cornfield 
and  the  wheatfield  turns  out  to  be  the  same;  that 
is,  it  tends  to  decrease  with  an  increase  of  reference 
wind.  U0.  The  corresponding  S  values  of  wheat 
are  about  10  times  those  of  corn.  The  reason 
for  this  difference  is  that  the  effective  area  density, 
.«4o,  of  the  wheatfield  is  much  higher  than  that  of 
the  cornfield. 

In  figure  16  the  measured-canopy  flow  profiles 
for  corn  and  for  wheat  are  plotted.  They  fit  well 
with  the  computed  wind  profiles.  This  validates 


Figi  re  !7. — Velocity  gardients  at  peak-plant  height  as 
a  function  of  shape  factor  S. 

the  assumption  of  the  turbulent  mixing-length 
concept  used  in  the  analytical  model. 

For  both  the  wheatfield  and  the  cornfield,  the 
ratio  of  eddy  viscosity  at  peak-plant  height  to 
the  viscosity  of  air  varies  from  1,000  to  8,000 
within  the  range  studied.  The  assumption  of  a 
fully  turbulent  model,  and,  consequently,  the 
application  of  mixing-length  theory  to  the  study 
ot  air  motion  in  the  vicinity  of  a  crop,  are,  there¬ 
fore,  valid. 

Shear -Wave  Theory 

During  the  course  of  this  study,  a  very  crude 
shear-wave  theory  was  proposed  to  study  one 
phase  of  the  nonsteady  state  phenomena.  This 
idea  was  generated  by  field  measurements  of  the 
wind  profiles  when  it  was  observed  that  wind 
components  at  different  heights  exhibited  phase 
lags,  thereby  suggesting  a  propagating  shear-wave. 
Figure  10  show's  tw'O  simultaneous  recordings  of 
the  horizontal  wind  at  two  vertical  positions — 
Z= 560  cm.  and  Z=300  cm  — within  a  cornfield. 
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The  lines  containing  high  frequencies  are  instan¬ 
taneous  recordings  from  anemometers  with  a 
response  time  of  less  than  one  second.  The 
second  set  of  lines  in  long  dashes  represents  the 
running  average  of  15  seconds,  and  the  third  set 
represents  a  running  average  of  60  seconds.  For 
the  60-second  running  average,  all  wavelengths 
under  60  seconds  are  averaged  out,  and  only  the 
wave  spectrum  greater  than  60  seconds  remains. 
These  long-period  waves  vary  nearly  in  phase 
with  each  other,  and  are  defined  as  quasi-steadv 
flow.  The  15-second  average  line,  which  oscillates 
about  the  60-second  average  line,  represents  the 
spectral  components  between  15-second  and  60- 
second  waves. 

These  waves  are  seldom  in  phase  between  the 
two  levels  under  observation.  The  top  wave 
may  either  lead  or  lag  behind  the  lower  waves. 
It  is  not  easv  to  differentiate  which  is  leading 
or  lagging  wfien  waves  of  multiple  periods  are 
intermixed.  However,  when  there  is  a  sudden 
increase  in  wind  intensity,  the  wind  on  the  upper 
level  generally  can  be  observed  to  be  leading 
that  of  the  lower  level.  This  type  of  oscillating 
wind  is  identified  as  nonsteadv  wind. 

The  differential  equation  for  shear-waves  prop¬ 
agating  in  the  Z  direction  can  be  obtained 
through  the  modified  Navier -Stokes  equation  by 
replacing  eddy  viscosity  for  molecular  viscosity 
and  neglecting  pressure  gradients  and  velocity 
gradients  in  the  i  and  y  directions. 


dt  dZ 


l>IKi 


Differentiating  equation  75  with  respect  to  Z 
and  multiplying  each  side  by  «(Z),  we  have 

d  W  t(Z)  d*r 

dtt{  ’dZ  p  dZ 1 
or 

|-K(Z)  ^  (76) 

The  eddy  diffusivity  K(Z)  may  be  taken  as 
linearly  proportional  to  Z;  however,  this  leads  to 
a  rather  complex  solution  for  the  shear  wave.  To 
simplify  the  situation,  it  is  expedient  to  consider 
that  K(Z)  is  a  constant,  with  an  average  value 
of  K0  within  the  range  considered.  Equation  76 
is  then  reduced  to 


dr  -.  d*r 


t=t04-t'  sin(»r). 

The  steady'  state  solution  for  equation  77  is 
similar  to  the  well-known,  periodic-heat-transfer 
problem  of  a  flat  slab  (2). 


^  ,  r'  cosh  w'Z(l  +  i) 

t(Z,I)  =  to+ - r — /v~  /.  i  -\  ’  (78) 

cosh  w  Z10(  l+i ) 


where 


■mat 


When  the  frequency  factor  w'Z10  is  less  than 
0.5,  the  absolute  magnitude  of  oscillating  shear, 
r'(Z),  differs  only  2  percent  at  Z=0  from  r'  at 
Z,0.  The  corresponding  phase  lag  at  these  levels 
is  14°.  That  is,  when  the  frequency  factor  is 
less  than  0.5,  constant-shear  condition  can  be 
assumed. 

It  can  be  shown  that  the  above  shear  wave  will 
propagate  downward  with  a  speed  of 


'  dl"  cosh  u'Zfl+i)  ”1  (79) 

dzLUrgcosh  W'Z10(1  +  ?)J 

The  speed  of  wave  propagation  is  not  constant 
with  respect  to  height.  However,  it  can  be 

shown  that  for  a  good  approximation,  V,  may  be 

taken  as 

Vr*=\2K^>,  ^>0.2.  (80) 

^10 

The  shear-wave  length  in  the  Z  direction  can 
be  defined  as 

Xr=^=2  (80 

w  M  w 

The  ratio  of  Z)0  to  X,  is 


-10  1  y  /  w  _Z,0O,' 

_ 2tt  2Kq  2ir  ‘ 


with  the  boundary  conditions: 


For  the  limiting  case  of  constant  shear, 
dr/dZw  0,  where  Z,0«'  is  less  than  0.5,  the  ratio  of 

^-°<?-=0.08. 

~  2 n 

That  is,  when  the  shear-wave  length  is  more  than 
12  times  the  reference  height,  Zl0,  constant-shear 
condition  exists  (quasi-steady  state).  The  shear- 
wave  length  required  for  this  state  of  flow  is 


and  r  at  reference  height  Z|fl  oscillates  about  a 
mean  shear  r0  as 


x’-aos=125ml',ere' 

From  equation  81  the  frequency  of  oscillation,/, 
for  quasi-steady  condition  can  be  obtained. 
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/  4  tK, 


(83) 


For  a  typical  condition  where  U0= 500  cm. /sec., 
A,  =  12,500  cm.,  and  K0=  50,000  cm.Vsec.,  the 
corresponding  value  of  /  from  equation  83  is 


,  4x50,000 
12,500* 


4X10  3  cyclcs/sec. 


Under  the  conditions  when  the  time  and  space 
coordinates  are  interchangeable,  the  equivalent 
velocity  wavelength  can  be  defined  as 


500 

4X10"3 


=  125,000  cm.  =  1,250  meters. 


Hence,  from  the  very  crude  shear -wave  theory,  an 
approximate  boundary  value  for  quasi-steady 
state  condition  is  defined.  That  is,  any  wave¬ 
length  greater  than  1,000  meters  may  he  considered 
as  a  quasi-steady  wind. 

The  wind  spectrum  can  now  be  divided  into 
the  following  classification  of  eddy  scales: 

Scale  Ratio  Eddy  Classification 

=^■<4. . Turbulent  eddies. 

^10 


4<^-<100_.  _  Nonsteadv  wind. 
^10 

100<  -=?•  - .  Quasi-steadv  wind. 

^10 


The  relaxation  time,  tr,  of  a  boundary  layer, 
referred  to  at  the  beginning  of  this  section,  is 
defined  as  the  time  for  the  shear  wave  to  propagate 
from  the  reference  height  down  to  the  plant. 


2,o — h 

VT 


(84) 


For  the  typical  example  as  given  above,  the  value 
for  tr  is  14  seconds.  Winds  of  higher  frequency  or 
shorter  wavelength  will,  of  course,  give  shorter 
periods.  The  phase-time  relationship  of  waves 
of  30  seconds  duration,  as  shown  in  figure  10  (in 
long  dashes),  is  about  4  seconds  between  the  two 
recording  levels.  Equation  84  gives  a  relaxation 
time  of  about  the  same  magnitude.  This  checks 
well  with  the  experimental  data. 


ON  THE  DIFFUSION  OF  HEAT  AND 
MASS  IN  A  TURBULENT  SHEAR 
FLOW 


By  D.  A.  Sp  KNCK 

The  feature  that  distinguishes  turbulence  from 
laminar  flow  is  its  randomness:  In  any  given 


case  it  is  possible  to  give,  at  the  very  best,  only  a 
statistical  description  of  the  way  m  which  the 
distributions  of  velocity,  pressure,  temperature, 
and  density  vary  about  their  means.  Fortu¬ 
nately,  such  a  description  is  usually  good  enough, 
since  one  is  only  interested  in  tne  behavior  of 
properties  averaged  over  a  time  that  is  long 
compared  with  that  of  a  turbulent  fluctuation. 
Exactly  as  in  statistical  mechanics,  one  is  only 
interested  in  the  mean  properties  of  an  ensemble, 
and  one  not  only  could  not  but  would  not  want 
to  say  how  an  individual  particle  behaves. 

The  analogy  with  statistical  mechanics  is  a 
natural  one,  and  led  Prandtl  to  introduce  the 
notion  of  a  mixing  length  as  a  counterpart  of  the 
mean  free  path,  but  it  is  nowadays  generally 
agreed  that  although  this  concept  has  dimensional 
validity,  it  is  a  great  oversimplification  of  the 
transfer  processes  at  work  in  a  turbulent  flow. 
For  this  reason,  one  is  on  dangerous  ground  in 
trying  to  predict  the  results  of  a  new  experi¬ 
ment — for  example,  in  applying  suction  or  blowing 
to  a  turbulent  boundary  layer — by  means  of 
mixing  length  theory,  even  though  it  can  be  used 
very  successfully  to  rationalize  some  well-known 
observations,  the  most  notable  being  the  loga¬ 
rithmic  law  of  the  wall. 

Turbulence  presents  a  much  more  difficult 
problem  than  the  mean  motion  of  molecules  in 
a  gas,  however,  since  the  motion  of  one  fluid 
particle  affects  that  of  every  other,  and  the 
simplification  of  “molecular  chaos”  is  not  avail¬ 
able.  Moreover,  the  governing  equations  are 
nonlinear,  and  the  fluctuating  motion  is  always 
three-dimensional,  even  in  a  one — or  two — di¬ 
mensional  mean  flow.  For  this  reason,  substan¬ 
tial  progress  in  predicting  the  development  of  a 
field  of  turbulence  from  a  statistically  known 
initial  state  has  only  been  possible  in  the  simplest 
case:  that  of  decaying  homogeneous  isotropic 
turbulence.  Definitive  accounts  of  this  subject 
are  given  by  Batchelor  (/)  and  Lin  (19),  and 
mention  should  also  be  made  of  the  review  by 
Liepmann  (16). 

The  primary  aim  of  workers  in  this  field  is  to 

(>redict  the  form  and  time  variation  of  the  turbu- 
ent  energy  spectrum.  In  this  way  they  obtain 
some  insight  into  the  mechanisms  by  which  energy 
is  transferred  by  the  inertia  terms  in  the  equations 
of  motion  from  large  eddies  with  low-wave 
numbers  through  an  intermediate  range  to  the 
high-wave  numbers  at  which  dissipation  takes 
place  by  viscous  action.  The  well-confirmed 
existence  of  an  inertial  subrange  in  which  the 
spectrum  depends  only  on  the  rate  of  dissipation 
and  the  wave  number,  k—  and  is,  therefore,  found 
by  a  simple  dimensional  argument  to  vary  as 
k~il3 — shows  that  the  process  is  a  continuous  one; 
i.e.,  that  energy  is  passed  successively  down  the 
whole  wave-number  scale  without  bypassing 
any  part  of  the  range.  To  make  further  prog- 
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ress,  however,  some  physieal  hypothesis — for 
example,  t hat  the  velocity  correlations  of  a  certain 
order  arc  those  given  by  joint  normal  distribu¬ 
tions — must  be  made,  and  nothing  completely 
satisfactory  seems  to  have  emerged,  so  the  theory 
of  homogeneous  turbulence  is,  at  present,  some¬ 
thing  of  a  standstill. 

In  the  meantime,  however,  t lie  same  methods 
have  l)een  used  to  aequire  a  considerable  insight 
into  the  qualitative  nature  of  turbulenee  in  the 
important  elass  of  flows  in  which  the  mean  motion 
is  sheared,  as  between  fixed  boundaries  in  a  pipe 
or  channel,  free  boundaries,  or  in  wakes  or  jets, 
and  one  fixed  and  one  free  in  a  boundary  layer. 
A  comprehensive  survey  of  these  flows  was  given 
by  Townsend  (37)  in  1956.  From  his  work  a 
fairly  definite  picture  has  emerged  of  the  way 
eddies  compare  in  size  with  the  width  of  the  layer 
extract  energy  from  the  mean  flow  and  hand,  it 
down  through  eddies  of  progressively  smaller  size 
to  those  in  which  it  is  ultimately  dissipated.  (The 
latter  are  not,  as  was  once  thought,  necessarily 
isotropic  even  in  a  shear  flow.)  The  net  rate  of 
production  of  turbulent  energy  per  unit  mass  is 
5r 

u  —i  the  product  of  mean  velocity  and  shear  stress 

gradient,  and  if  it  were  possible  by  considering 
the  structure  of  the  turbulenee  to  say  how  this 
quantity  depends  on  the  local  velocity  and  pres¬ 
sure  gradients,  the  calculation  of  mean  flow  as  a 
shear  layer  would  be  fully  determinate;  a  further 
relation  being  provided  by  t  lie  Navier-Stokes  equa¬ 
tions.  But  Townsend  and  his  coworkers  are  still 
quite  a  long  wav  from  being  able  to  make  any 
(juantitative  predictions  of  the  transport  properties 
in  terms  of  local  mean  quantities.  (“Transport” 
is  used  so  as  to  include  heat  and  mass  diffusion 
into  shear  force;  i.e. ,  momentum  transport,  in  the 
discussion.) 

Actually,  it  is  not  at  all  clear  that  a  representa¬ 
tion  in  terms  of  local  quantities  is  physically 
meaningful,  even  though  it  is  always  formally  pos¬ 
sible  by  means  of  an  eddy  diffusion  coefficient. 
Transport  of  any  quantity  is  brought  about  to 
different  extents  by  eddies  of  different  sizes.  In¬ 
tuitively  one  has  a  picture  of  a  lump  of  fluid 
carrying  with  it  some  property  having  the  mean 
value  appropriate  to  its  initial  surroundings  and 
being  eonveeted  by  an  eddy  to  a  point  where  the 
mean  properties  are  different,  and  where  it  then 
shares  the  transported  property  with  t lie  surround¬ 
ing  fluid.  The  distance  traveled  by  the  lump  is 
then  the  mixing  length,  and  the  flux  of  the  trans¬ 
portable  property  should  bo  proportional  on  this 
argument  to  the  mean  gradient  over  the  length. 

But  lumps  converted  by  different  eddies  will 
travel  different  distances.  At  least  in  the  case  of 
the  layer  eddies,  an  exchange  coefficient  based  on 
the  local  gradient  at  one  point  in  the  shear  layer 
does  not  seem  appropriate,  both  l>ocause  the  eddies 
are  comparable  in  width  to  the  shear  layer  and 


because  the  exchange  process  takes  place  con¬ 
tinuously  by  molecular  action  wherever  there  is 
a  local  concentration  gradient,  not  waiting  for  the 
lump  to  reach  its  final  position.  In  fact,  Town¬ 
send’s  experiments  (see  Hinze  (.9),  pp.  288-289) 
indicate  that  momentum  transport  is  primarily 
due  to  gradient-type  diffusion  by  small-scale  tur¬ 
bulence  (presumably  because  the  u  and  v  compo¬ 
nents  in  the  large  eddies  tend  to  be  uncorrelated), 
whereas  heat  is  transported  both  by  gradient  dif¬ 
fusion  and  by  the  bulk  mechanisms  associated 
with  large  eddies. 

Thus,  as  pointed  out  by  Lighthill  (18)  in  review¬ 
ing  Townsend’s  book,  it  is  physically  quite  signifi¬ 
cant  to  define  an  eddy  viscosity,  e,  Dv 


—pit'v'  =  ( 


(yu 


and  to  expect  t  to  depend  on  local  quantities.  To 
this  is  added  the  viscous  contribution  to  shear 
stress,  giving  altogether 


t=(^+ *) 


du 

dy 


and  the  analogy  between  the  gradient  diffusion  due 
to  eddy  and  molecular  motions  is  clear.  By 
analogy  with  the  first  equation,  an  eddy  con¬ 
ductivity  K  and  an  eddy  diffusion  coefficient  D 
(for  a  primary  mixture)  may  be  defined  in  terms 
of  the  fluxes  of  enthalpy,  h,  and  mass,  denoted  by 
concentration,  c,  of  one  species,  due  to  the  turbu¬ 
lent  motion  by 


if  ,  dA  —7—7  y\  dc 

— ph  v  —K  ■=-->  — pc  r  —D  5- — 

t>y 


The  corresponding  molecular  transport  coeffi¬ 
cients  are  denoted  by  k  and  D,  say.  (Note  that 
the  specific  heat  C„  has  been  absorbed  in  the 
definitions  of  k  and  K.)  The  kinetic  theory  for  a 
dilute  gas  at  room  temperatures  shows  that  the 
ratios 


are  all  of  order  unity  (between  0.7  and  1 .4,  roughly) 
and  one  is,  therefore,  able  to  treat  the  diffusion  of 
heat  or  mass  by  analogy  with  that  of  momentum. 
For  example,  Reynold’s  analogy,  wdiich  holds 
exactly  for  7V=1,  states  that,  when  suitably  non- 
dimensionalized,  skin  friction  is  equal  to  heat 
transfer  for  a  boundary  layer. 

For  turbulent  flows  one  naturally  seeks  to  use 
the  same  type  of  analogy,  even  though  for  the 
reasons  discussed,  it  rests  on  less  certain  physieal 
grounds.  Nevertheless,  it  is  possible  to  give  quite 
good  quantitative  answers  to  practical  problems 
in  chemical  engineering,  micro-meteorology,  etc., 
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by  relating  the  fluxes  of  teat  and  matter  to 
the  shear  stress  by  means  of  the  corresponding 
parameters 


(Pr) 


which,  in  view  of  the  broadly  analogous  exchange 
mechanisms,  should  all  be  somewhere  near  unity. 
General  surveys  of  this  approach  have  been  given 
by  Ijees  (IS)  and  by  Spalding  (SO),  the  latter  being 
particularly  concerned  with  combustion  problems. 
The  usual  practice  is  to  give  the  ratios  constant 
values  in  any  given  shear  flow.  In  this  way  the 
present  author  (SI)  was  able  to  relate  the  distribu¬ 
tion  of  enthalpy  to  that  of  shear  stress,  both  being 
regarded  as  functions  of  velocity,  in  terms  of  an 
arbitrary  value  a  for  the  turbulent  Prandtl  num¬ 
ber.  If  the  energy  equation  was  integrated,  it  was 
then  possible  to  express  both  the  recovery  and 
Reynold’s  analogy-factors  in  terms  of  a,  and  by 
talang  of=?0.85,  which  is  midway  between  unity 
and  ( Pr )  moucui,r,  good  agreement  with  the  known 
experimental  results  for  air  was  found. 

In  view  of  the  differences  in  transport  mech¬ 
anisms  for  heat  and  momentum,  it  would  be 
surprising  if  Prmrb  actually  were  constant  across 
the  layer  and  the  above  value  should  probably 
be  looked  on  as  a  weighted  means.  The  transport 
must  depend  on  the  eddv  structure,  and  this 
varies  across  the  layer.  For  the  region  where 
small  eddies  predominate  (close  to  the  sublayer), 
gradient  diffusim  would  probably  be  the  main 
mechanism,  and  one  might  expect  to  find  ( Pr 
~(Pr)moUeuUr,  but  near  the  outer  edge  of  a  shear 
layer  large  eddies  are  predominant  and  transport 
is  primarily  by  bulk  motions,  so  it  seems  more 
probable  that  the  Prandtl  number  should  then 
approach  unity.  The  integration  of  reference 
(SI)  can  actually  be  carried  out  formally  for  an 
arbitrary  variation  of  ( Pr ^)tMrtl,  as  is  shown  in  the 
notation  below,  giving  a  result  of  some  generalitv 
for  zero-pressure  gradient  flows.  Like  most  such 
results,  they  are  not  suitable  for  application 
without  some  physical  hypothesis.  In  this 
case,  the  dependence  of  (Pr)tVT-,  on  the  shear  stress 
may  be  a  suitable  representation  of  the  local  eddy 
structure. 

The  same  methods  can  be  used  to  estimate  the 
diffusion  of  mass  across  a  shear  flow  in  terms  of 
shear  stress  and  velocity  in  terms  of  an  arbitrary 
Schmidt  number,  and,  in  the  absence  of  definite 
experimental  information,  one  would  take  a  value 
somewhere  between  unity  and  the  molecular 
value  as  the  most  likelv.  Meteorologists  usually 
use  mixing  lengths  having  some  prescribed- 
height  dependence  to  obtain  the  distribution  of 
velocity  and  shear  stress,  as  detailed  bv  Priestlev 
(25).  When  this  is  known,  one  can  proceed 
exactly  as  in  the  beat  diffusion  case  to  calculate 


the  flux  of  a  foreign  gas  throt  gh  the  atmosphere 
due  to  turbulent  mixing. 

Integration  of  the  energy  equation  in  a  zero- 
pressure  gradient,  with  variable  turbulent  Prandtl 
number. — The  energy  equation  with  the  usual 
notation  is 


dh  .  .  ,  —ri\  dh  dq  , 

pU_+(pP+pt,  , 


du 

by 


and  if  h—h(u),  this  can  be  written, 


1  dr  dh _ 1  rfg  j 

r  du  du  r  du' 

The  heat  flux, 


«=<*+*>  !=*••$ 


say,  and  likewise 


ini 

T=™by 


If  we  write  ^ =a ,  then  *  is  a  function  of  the 

«ell 

Irtarb~k~kttt-k 
In  terms  of  a  the  energy-  equation  is 


d*h 
du 1 


+[('-»)  ±  In  (£)-^  (ln«)]^— «• 


(The  wall  stress  tk  has  been  brought  in  for  con¬ 
venience.) 

Denote  by  <j>  the  Stieltjes  integral 

Then,  if  we  set  a—o  at  the  wall  u=0,  the 
result  of  integrating  the  energy  equation  twice  is 


If  a = constant,  this  reduces  to  the  usual  Crocco 
integral. 

It  can  also  be  integrated  in  a  fairly  straight¬ 
forward  manner  if  a  is  assumed  to  be  a  simple 

function  of  7  ■  Forinstance.it  might  be  physically 

T  it 

realistic  to  assume  a  linear  dependence,  by  setting 


a=  1  —  ( 1  —  a ) 


which  goes  to  unity  at  the  free  boundary  and  to 
a  at  the  wall. 


A  METHOD  FOR  THE  COMPUTATION  OF  LOGARITHMIC  WIND 
PROFILE  PARAMETERS  AND  THEIR  STANDARD  ERRORS 


By  Wtnton  Covey 


The  logarithmic  wind  profile  equation  may  be 
written: 

r.v.i„(£±»>  <85) 

Here,  z  m  height  above  the  ground,  Fis  mean  wind- 
speed  ov  r  a  short  (e.g.,  15-minute)  time-interval. 
Equation  85  is  observed  to  fit  observations  well 
when: 

(1)  The  anemometry  is  good 

(2)  Thermal  stratification  is  adiabatic, 

(3)  The  measuring  site  is  amid  a  uniform 

patch  of  ground, 

(4)  The  lowest  anemometer  is  not  too  low, 

and 

(5)  The  highest  anemometer  is  not  too  high. 

Requirements  3  and  4  lead  some  people,  in¬ 
cluding  the  author,  to  view  the  logarithmic  profile 
as  applying  basically  to  area-average  profiles, 
rather  than  to  point  profiles.  Requirements  3 
and  4  are  then  looked  upon  as  requirements 
for  obtaining  area-average  profiles  with  a  minimum 
of  instrumentation. 

The  interpretation  put  upon  the  three  param¬ 
eter? — F»,  z0  (sometimes  called  roughness  length), 
and  D  (sometimes  called  displacement  height) — 
has  been  extremely  dn  >rse  during  the  past  24 

} rears.  Extreme  viewk  ave  been  that  roughness 
ength  is  "a  constant  of  integration”  on  the  one 
hand,  and  a  physical  parameter  of  the  ground 
surface  on  the  other.  The  parameter  F*  has 
generally  been  identified  as: 


Recent  work  shows  that  for  flexible  elastic  plant 
covers  such  as  corn  or  rice,  both  z0  and  D  vary 
with  windspeed,  as  indicated  by  F*,  geostrophic 
wind,  or  speed  at  a  reference  height.  The  study 
of  the  manner  of  this  dependence  would  be  facili¬ 
tated  by  a  convenient  means  of  determining  F*, 
Zo,  D,  and  their  standard  errors  from  observations. 

On  receipt  of  some  notes  from  Stephen  M. 
Robinson  and  C.  B.  Tanner  4  at  the  University 
of  Wisconsin,  giving  substantially  the  material 
in  pages  1  to  5  of  Technical  Note  No.  1  {26), 
a  program  for  machine  computation  of  adiabatic 
wind  profile  parameters  was  written  here.  This 
program  accomplishes  the  same  data  reduction 
task  as  Mr.  Robinson’s,  and  is  based  on  equations 
written  bv  him.  It  has  two  major  differences 
in  technique  that  may  prove  useful.  The 
subprogram  for  standard  errors  and  the  equa¬ 
tions  on  which  they  are  based  are  original  con¬ 
tributions.  The  program  is  written  in  the  Bur¬ 
roughs  version  of  Algol,  for  the  Burroughs  220 
computer.  This  report  discusses  mathematical 
methods  and  not  the  details  of  the  computer 
program. 

REVIEW  OF  THE  PROBLEM 

The  equation  to  be  fitted  to  the  adiabatic 
wind  profile  observations  is: 

F=F*  In  (““}  (85) 


(86) 


where  k  is  Von  Karman’s  constant,  ~0.4;  r 
is  shear  stress  in  the  surface  layer;  and  p  is  air 
density.  Let  z,  be  the  value  of  z  at  which  equa¬ 
tion  85  gives  F=0.  Then: 


or, 


0=F*  In 


z,+D 

- -  f 

*0 


z,  =  —  I)-\-Z(,.  (87) 


The  parameter,  I),  is  therefore  related  to  the 
height  at  which  the  profile  extrapolates  to  zero, 
differing  by  c0,  and  having  a  change  in  sign. 

For  rough,  rigid  surfaces,  z0  appears  to  be  a 
physical  parameter  of  the  surface  and  z,  =  0. 
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Robinson’s  development  is  followel  in  general, 
with  summation  and  mean  value  notation  rather 
than  vector  notation.  Data  are  speeds  V,  for 
the  n  measuring  heights,  zt. 

The  sum  of  the  squares  of  errors  to  be  minimized 
in  evaluating  F, ,  D,  and  z0  is: 


i^et : 


W=ll)  z0, 
j(=ln  ( zt-f-D ), 


1  " 


(881 


4  Personal  communication. 
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yt=xt-x, 

ut=Vt-V, 

rt=\l(2f+D). 

Put  the  three  partial  derivatives  simultaneously 
equal  to  zero : 


— 2£  qr(-r,(/,-tr)K*,-w)}=o;  m 

ov  *  tTj 

{ir,-X’,(x1-«;)]r<}=0;  (90) 

OU  1-1 


A  form  of  the  latter  function  is  chosen  by 
Robinson.  The  choice, 


?sF  IV  n 

^=£i*v{r(-p,(/(-w)}=o. 

02 o  20  f-1 


In  equations  89,  90,  and  91,  w  is  a  function  of  z0 
and  the  x(  and  r<  depend  upon  D.  The  set  of 

values  (I7*,  D,  z0)  is  required  that  satisfies 

equations  89,  90,  and  91. 

In  equation  91,  averages  are  taken,  and  the 
equation  is  rearrangeu  to  give: 

w=x-y~-  (92) 

Equations  90  and  89  are  rewritten  and  averaged, 
and  the  value  of  w  from  equation  92  is  substituted 
in  them.  This  gives: 

0=ur— U»  ~yr  (93) 

and 

0=u7/— fTy3.  (94) 

From  equation  93  comes  the  equation  for  U*: 

u,==-  (95) 

Combine  equations  94  and  95  to  obtain: 

0=uT/-^5-  (96) 

yr 

Equation  96  is  the  implicit  equation  in  the  single 
unknown,  D,  which  is  to  be  solved.  This  is  done  by 
defining  a  function,  g(D),  which  has  as  a  single 
real  root  the  desired  value  of  D,  and  by  finding 
that  root.  Any  of  several  functions  g(D)  will  do: 

y'(/?)=constant  •  {uy— p  ur/ijr } 
or 

f/"(/>)=con8tant  •  { y*  ur  —  uij  1/r}, 
for  example. 


<KZ?)=constant 


14-^4 

\  </>  uy-U'J 


has  the  merit  that  g(D)  so  defined  has,  for  reason¬ 
able  data,  a  nearly  linear  form  for  a  wide  range  of 
trial  values  of  D.  This  near-linearity  is  convenient 
when  new  estimates  of  the  root,  D,  are  found  bv 
linear  interpolations.  Hence,  equation  97  is  used. 

The  second  difference  in  technique  is  in  finding 
the  first  two  and  subsequent  estimates  of  D.  In¬ 
stead  of  a  pair  of  reasonable  guesses  being  supplied, 
the  program  automatically  makes  two  extreme 
guesses,  to  bracket  the  root.  New  estimates  are 
found  alternately  by  linear  interpolation  and  by 
bisection,  until  one  of  three  criteria  calls  for  making 
a  final  linear  interpolation.  Each  new  estimate  is 
used  to  narrow  the  bracket. 

The  data  format  consists  of  integers  (groups  of 
digits)  separated  by  spaces: 

aaa  bb  ccc  ccc  ccc  ccc. 

bb  integers 

If  the  first  integer  (aaa)  ^ 0 ,  then  (aaa)  is  the 
profile  identification  number,  and  the  (ccc)  are 
velocities,  starting  at  the  lowest  anemometer 
height.  The  second  integer  (bb)  is  the  number  of 
measuring  heights.  If  (aaa)=0,  then  the  (ccc) 
are  anemometer  heights,  beginning  with  the  lowest. 
The  computer  subsequently  converts  input  data  to 
floating  point  form,  and  the  output  has  explicit 
decimal  points.  Units  of  I7,  in  the  output  are 
those  of  V,  in  input;  units  of  Z?and  z0  in  the  output 
are  those  of  zt  in  input.  Units  of  centimeters  per 
second  and  centimeters  were  visualized  when 
writing  the  program,  but  are  not  essential. 

The  first  and  second  trial  values  of  D  are 
D]=  —Z]-\-l  and  J?,=  +  Z|  —  1.  If  these  do  not 
establish  a  bracket  for  D,  the  fact  is  noted  and  the 
next  profile  is  read  in. 

The  three  criteria  for  making  a  final  estimate  of 
D  by  linear  interpolation  are: 

(1)  The  absolute  value  of  g(D)  less  than  a 

critical  amount,  or 

(2)  The  size  of  the  bracket  less  than  a  critical 

amount,  or 

(3)  The  number  of  trial  values  of  D  reaches  a 

maximum  allowed  number. 

With  several  sets  of  real  data,  all  profiles  having 
four  or  more  measuring  heights,  and  heights  being 
measured  from  approximately  the  mean  soil  sur¬ 
face,  seven  trial  values  of  D  (including  the  first 
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two)  were  always  sufficient  to  evaluate  the  d»-  t(em.)  V  (mm./sec.) 

placement  height,  D,  to  the  nearest  hundredth  of  40 _  .4, 902 

a  centimeter,  if  the  least-squares  value  was  within  80. _ . . . _ _  5,  853 

the  original  bracket  (fig.  18).  160 _ _ _ 6,737 

320 _ _ _ 7,592 

A  SAMPLE  COMPUTATION  540 -  8,375 


.  Successive  estimates  of  D  and  corresponding 
The  same  data  are  used  as  in  Robinson  s  sample  values  of — 

calculations.  These  are  Lettau’s  mean  adiabatic  _ 

profile  from  the  Johns  Hopkins  observations  at  n(P)=L—  — 

O'Neill,  Nebr.,  in  1953  (15).  if  tiyyr 


Method:  D  9(D)  Resulting  bracket 

-39.  0  +  0.  029518  . 

+  39.0  -  0.060772  -  39.0  +  39.0 

By  linear  interpolation... . .  -13.499818  +  0.006182  -13.5  +  39.0 

By  bisection _ _  +12.750091  -0.030611  -13.5  +12.75 

By  livear  interpolation . . —9.088715  —  0.000678  —13.5  —  9.09 

By  bisection . -11.294266  +  0.002713  -11.29  -  9.09 

By  linear  interpolation . -9.529773  -  0.000006  -11.29  -  9.529773 


Then  a  final  linear  interpolation  gave  —9.53  cm.  approximately  linear  dependence  of  g(D)  on  D  for 

for  D.  ^"hc  corresponding  2#  was  0.432  cm.,  and  the  example  above. 

V*  was  1,151  mm./sec.  Figure  18  shows  the 


ESTIMATING  THE  ACCURACY  OF  F*, 
D,  AND  z0  FROM  ADIABATIC  WIND 
PROFILE  DATA 


dD 

av,  v. 


dH\ 

aVi  v|t  v, - p 

dH 

dD  v„  r. 


(100) 


Using  a  computer  with  eig.,t  significant-figure 
arithmetic  does  not  lead  to  computational  error 
in  the  profile  parameters.  Errors  depend  upon 
the  data  for  the  heights,  zu  and  the  s  peeds,  V,.  The 
heights  can  be  assumed  to  be  accurate,  and  all  the 
error  of  measurement  (sampling  and  instrumental) 
can  be  assigned  to  the  speeds.  Then  statistical 
analysis  leads  to  estimates  of  the  errors  in  com¬ 
puted  parameters  V*,  D,  and  ^  as  functions  of 
the  data.  This  statistical  analysis  involves  only 
the  assumptions  already  implied  in  the  fitting  of  a 
least-squares  profile. 


Define 


The  Error  in  D 


H(D,  Vu  Vj,  .  .  .  Vn)=uy  yr-if*  ur  (98) 

where  D  is  an  arbitrary  variable,  and  the  V,  are 
measured  velocities. 

Then, 

Divide  equation  99  by  dVu  and  put  0=dVj= 
dV,=  .  .  .  dV,=dH.  Then: 


and,  of  course,  similar  equations  apply  with  respect 
to  V,.  \\  .  .  .  Vn.  The  error  in  displacement 
height,  D,  is  given  bv: 

Em~w.  '':<v')+m  E(Vt)+  ■  +sv.  E(V-> 

+  ( higher  order  terms) 

where  the  derivatives  are  evaluated  at  the  true 
values  of  the  velocities. 

The  standard  error  of  D  is  given  (approximately) 

'TO={g(||Wj}W  001) 

If  the  standard  error  of  VT<  is  assumed  constant 
with  height,  then: 

»«»-'-!?  {s (S’}"  «“> 


Equations  are  required  for 


<  |-v  ttHQ  p.  rr 

v,.vt _ r„  )  v,*r  d 
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Therefore, 


dF»  yr—vr  rs+ur  y*  dD  + 

*>V,  y,±,~  (' W )*  dVjV,^+nyr 

(107) 

Evaluating  at  the  true  value  of  velocities, 


;  _ _ t 7  •  ° 

at/  v  *  rrr 


ra  0D 


yrdVjIv,^,  nyr 

The  Error  in  z0 

The  roughness  length,  z0l  is  given  by: 

z0=exp(x—~y  (108) 

The  standard  error  of  z0  is  approximated  by 


SE(z0)~ 

But 

• 

dz0 

dFj  v,*, 

This  leads  to 

dz0 

Si 

dF, 

dz0 1  dD  I 


«  NT? 


?+(fM-ur*  yr-urr* 


+"r^)^L+F.[v'sF-;]}-  <H1> 

If  we  evaluate  at  true  values, 


r 

(.  V  Ts) 

dD 

*  V 

V  \\y r) 

1  dFj 

1  (tj  1\^ 

'  y  yv  nV*yr  n)  j  For  a  SPt  °f  good  adiabatic  wind  profdes,  six 

J  measured  at  O’Neill,  Nebr.,  in  the  summer  of 

To  compute  these  standard  errors  of  the  profile  1956  by  Halstead  and  associates  (6)  were  chosen, 

parameters,  several  intermediate  variables  are  Data  and  results  of  analysis  are  tabulated  below. 


required  in  addition  to  those  found  in  computing 
the  profile  parameters.  These  are,  for  the  least 
squares  values  of  D,  F*,  and  20:  r  Jr,y»,T*;  and. 
at  each  measuring  height,  i:  #,  =r,  —  r; 


5 D 

dF, 


from  equation  100; 


from  equation  107; 


from  equation  111. 


5F* 

6F, 


dZo 

dV,  v,„, 


The  standard  error  of  the  windspeeds  is  afoo 
required,  SE(V,).  This  may  be  estimated,  for 
ti >4,  by 

=V^=r  (112) 

In  any  case,  SE(Vf)  is  at  least  as  large  as 
known  anemometer  errors.  The  parameters  for 
the  sample  problem  are  now  repeated,  with  their 
standard  errors  indicated: 

\\  =  1 15.1  ±  1 .4  cm. /sec. 

20=0.4d±0.04  cm. 

D=  — 9.5±  1.2  cm. 

ANALYSIS  OF  A  SET  OF  NEAR- 
ADIABATIC  PROFILES 


Date 
Time 
Tum  —  T.  jjm 

U||m 

U»m 
U»m 
U|m 
U.;m 
U»m 
V*,  cin/sec 
z o>  cm. 
I),  cm 


10  July  11  July  23  July  24  July  24  July  24  July 

1905  0005  1905  0605  0705  1905 

+  0.24°C  +0.23  -0.07  +0  48  -0.92  -0.18 

729  cm/sec  780  cm/sec  731  cm/sec  584  cm/sec  843  cm/sec  503  cm/sec 
639  689  635  532  763  427 

569  607  568  505  708  398 

497  543  _  420  609  339 

432  474  437  380  552  291 

373  411  377  336  501  259 

299  340  310  289  426  224 

108  +  4  113  +  5  108  +  6  75  +  6  107  +  6  Not  computed 

2.  0  +  0.  5  1.7+0.  4  2.  0  +  0.  7  .  65±.  38  .61 +.26 

9. 0  +  5. 3  12.0  +  6.0  12.4  +  8.5  5.3+10.3  9. 6  +  8.  0 
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Compared  to  the  mean  profile  first  discussed, 
these  six  have  appreciably  higher  standard  errors 
for  the  parameters.  That  D-Zo  is  positive  is 
a  minor  puzzle.  It  indicates  that  the  windspeeds 
at  the  highest  levels  are  less  than  expected,  while 
those  at  the  lowest  levels  are  greater.  At  the 
highest  level,  the  log-law  for  area  averages  does 
not  hold  so  well;  at  the  lowest  levels,  the  measure¬ 
ments  are  not  sufficiently  representative  of  area 
averages.  In  this  instance,  the  tendency  to  avoid 
local  obstacles  and  depressions  in  locating  the 
anemometer  mast  appears  to  have  led  to  selecting 
a  spot  with  slightly  higher  winds  than  average 
near  the  surface.  (These  data  were  not  taken 
with  this  manner  of  statistical  profile  analysis  in 
mind.  Subjective  analysis  of  profile  data  gen¬ 
erally  weights  the  anemometer  readings  according 
to  expected  reliability.)  These  new  computations 
show  again  that: 

(a)  Site  selection  and  anemometry  were  good ; 


(b)  The  magnitude  of  displacement  height, 
D,  was  small  compared  to  median  height 
of  measurement;  and 

(c)  2o  was  approximately  1  cm. 

SUMMARY 

Another  program,  based  on  Robinson’s  equa¬ 
tions  (26),  has  been  written  for  evaluating 
logarithmic  wind  profile  parameters.  Two  fea¬ 
tures  in  the  technique  differ  from  Robinson’s:  (1) 
the  method  of  choosing  the  first  estimates  of  D, 
and  (2)  the  method  of  finding  subsequent  estimates 
of  D.  These  perhaps  are  improvements.  In 
addition,  this  program  computes  good  estimates 
of  the  standard  errors  of  the  three  parameters, 
provided — 

(a)  The  accuracy  of  the  anemometers  is 
known,  or 

(b)  A  sufficient  number  of  measuring  heights 
are  used. 


TURBULENT  TRANSFER  CHARACTERISTICS  OF  THE  AIRSTREAM 
IN  AND  ABOVE  THE  VEGETATIVE  CANOPIES  AT  THE  EARTH’S 
SURFACE 


By  Jerry  Stoller  and  Edgar  R.  Lemon 


Turbulent  exchange  calculations  for  sensible 
heat,  latent  heat  (evaporation  and  condensation), 
and  carbon  dioxide  (photosynthesis  and  respira¬ 
tion)  are  important  m  estimating  some  of  the 
components  of  the  energy  balance  near  the  earth’s 
surface.  Basic  to  these  calculations  is  the  deter¬ 
mination  of  the  transfer  coefficient  for  momentum 
(Km)  and  applying  it  to  the  transfer  of  other 
molecular  properties.  In  essence,  one  implicitly 
determines  the  value  of  Km  in  the  turbulent 
boundary  layer  above  the  surface  and  uses  it 
to  calculate  the  exchange  of  other  molecular 
properties,  knowing  their  potential  drop  (gradient) 
across  some  known  distance  in  the  same  boundary 
layer.  The  assumption  is  made  that  the  momen¬ 
tum  transfer  coefficient  is  quantitatively  equal 
to  the  other  transfer  coefficients,  i.e.,  for  heat, 
water  vapor,  and  carbon  dioxide. 

Upon  the  accurate  measurement  of  the  momen¬ 
tum  transfer  characteristics  of  the  turbulent 
boundary  layer,  then,  rests  the  whole  aerodynamic 
approach  to  the  energy  budget. 

TURBULENT  EXCHANGE  ABOVE  THE 
CROP  CANOPY 

Where  vegetation  exists  at  the  earth’s  surface, 
the  turbulent  boundary  layer  has  to  be  sub¬ 
divided  into  two  regions:  One  above  the  plant 
surfaces  where  no  sources  or  sinks  occur  and  the 
other  within  the  vegetative  canopy  where  sinks 
and  sources  do  occur.  Exchange  calculations 
above  the  canopy  rely  upon  the  log  wind  profile 
law  and  the  logarithmic  law  of  vertical  distribu¬ 
tion  of  the  other  physical  quantities.  Three 
specific  problems  arise,  however,  in  the  use  of 
this  aerodynamic  approach  to  exchange  calcula¬ 
tions  above  the  vegetative  canopy: 

(1)  The  errors  created  by  appreciable  thermal 
gradients  in  the  turbulent  boundary  layer; 

(2)  Instrument  and  interpretative  limitations 
at  low  mean  wind  velocities;  and 

(3)  The  functional  relationships  of  the  log 
wind  profile  characteristics  to  both  mean  wind- 
speed  and  surface  properties  (geometric  and  elas¬ 
tic  properties  of  the  roughness  element). 

When  the  temperature  of  the  air  near  the 
vegetation  rises  significantly  during  the  daytime, 
then  thermal  convective  transfer  processes  (“free 
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convection”)  become  important  in  addition  to  the 
mechanical  convective  transfer  processes  (“forced 
convection”).  [The  forced  convection  processes 
are  responsible  for  the  log  profile  relationship.) 
Free  convection  superimposed  on  forced  convec¬ 
tion  causes  an  underestimation  of  the  transfer 
processes  when  the  log  profile  law  is  used,  and  is 
especially  pronounced  at  low  wind  velocities  when 
the  forced  convection  component  is  small.  Some¬ 
times  the  mean  wind  velocity  is  reduced  to  such 
a  low  level  that  it  is  not  measurable  bv  standard 
means  (rotating  cup  anemometers).  Under  these 
conditions,  the  logarithmic  law  of  vertical  distri¬ 
bution  of  physical  quantities  no  longer  holds 
true  and  the  calculated  fluxes  erroneously  go  to 
zero.  The  subject  of  “free  convection”  error 
wdll  be  dealt  with  in  detail  in  a  later  report. 

In  addition  to  free  convection  and  low  wind- 
speed  problems  associated  with  the  log  law  ap¬ 
plication  to  transfer  processes  in  the  turbulent 
airstream,  another  problem  arises  out  of  the  func¬ 
tional  relationship  between  windspeed  and  the 
wind  profile  characteristics  over  vegetative  sur¬ 
faces.  This  problem  is  created  by  the  changing 
geometric  properties  of  the  surface  (plant  growth) 
and  the  changing  elastic  properties  of  the  surface 
(waving  of  stalks  and  leaf  flutter).  This  func¬ 
tional  relationship  will  be  taken  up  in  detail 
first. 

The  eddy  flux  across  a  plane  of  some  entity 
(i.e.,  momentum)  which  is  neither  destroyed  nor 
created  in  its  transfer  (not  strictly  true  for  momen¬ 
tum)  can  be  written: 

Q=K^,  (113) 

where  Q  is  the  total  amount  of  the  entity  passing 
downward  through  a  unit  area  per  unit  time;  K  is 
the  transfer  coefficient,  and  dc/dz  is  the  gradient 
across  the  plane.  The  amount  of  the  entity,  c, 
is  expressed  on  a  per-unit-volume  basis,  and  z  is 
the  axial  distance  normal  to  the  plane. 

For  the  flux  of  momentum  (shearing  stress), 
t,  the  basic  equation  becomes: 

T-=Km~p  (114) 

where  u  is  the  windspeed  and  p  is  the  density  of 
the  air. 
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In  the  turbulent  boundary  layer  immediately 
above  the  vegetation  surface,  the  vertical  gradient 
of  wind  is  expressed  by : 


/roVl 
dz~k\P)  z 


(115) 


and  upon  integration  gives  the  “logarithmic  wind 
profile  law’’  where  the  windspeed,  u:,  at  height 
z  above  the  ground  is  represented  by: 


l'1  wherez>20,  016) 


t0  is  the  shearing  stress  at  the  surface,  r0  is  the 
roughness  length,  and  k  is  the  Von  Karman 
constant  (0.4). 

Equations  1 15  and  116  imply  that  the  windspeed 
profiles  are  functions  of  log  height,  r,  above  a 
given  reference  plane.  As  the  ground  vegetation 
growB,  however,  the  turbulent  boundary  layer  is 
displaced  upward,  causing  the  reference  plane 
from  which  z  is  measured  to  be  displaced  upward 
also  from  the  ground  surface.  It  then  becomes 
necessary  to  introduce  a  zero  point  displacement, 
d,  and  rewrite  equation  116,  thus: 


Zi+D 

z0 


(117) 


where  D=d+z0,  or  the  “effective  displacement 
parameter,”  and  zt  is  the  nominal  height  above 
the  ground  surface. 

The  friction  velocity,  V*,  defined  by: 


then: 


and, 


Km=(z,+D)kV*. 


(118) 

(119) 

(120) 


The  preceding  section  (p.  28)  takes  up  the 
method  used  here  for  the  determination  of  the 
logarithmic  wind  profile  parameters,  V’*,  z0,  and 
D  from  sufficient  experimental  data  u  on  z,. 

In  that  section,  T 


Alfalfa  Investigations 

Numerous  wind  profile  measurements  over 
alfalfa  were  taken  with  Sonoya  Rotating  Cup 
Anemometers  on  July  29,  30,  and  31,  1960,  when 
the  alfalfa  was  75  cm.  high.  The  equipment  and 
site  have  been  fully  described  elsewhere  (14)- 
Mean  wind  velocities  were  measured  over  5- 
minute  periods  at  90,  100,  120,  and  160  cm.  from 
the  ground.  Sampling  times  were  selected  when 
it  was  felt  that  thermal  effects  upon  the  profiles 
would  be  negligible  (i.e.,  when  temperature 


gradients  were  small).  Three  conditions  mini¬ 
mizing  this  error  were  observed :  (1)  The  measure¬ 
ments  were  made  near  the  surface  (within  160  cm. 
of  the  ground,  or  85  cm.  of  the  crop  top);  (2)  the 
measurements  were  made  during  cloudy  periods 
when  the  radiation  load  was  small;  and  (3)  the 
alfalfa  was  lush  with  plenty  of  soil  moisture 
(soil  moisture  tension  was  0.2  atm  at  6-inch 
depth). 

Table  1  and  figure  19  present  the  functional 
relationships  between  the  mean  windspeed  at 
160  cm.  and  the  profile  parameters  z0,  D,  and  V*. 
It  can  be  seen  that  the  roughness  length,  z0.  ranged 
from  28.5  to  1.3  cm.,  and  the  effective  displace¬ 
ment,  D,  from  —18.8  to  —65.0  cm.  The  friction 
velocity,  V*,  appeared  fairly  independent  of  wind- 
speed.  (Note  that  table  1  and  figure  19  give 
values  of  V*/'k  instead  of  V*.)  This  would  indi¬ 
cate  that  the  shearing  stress  was  fairly  constant 
and  independent  of  the  windspeed  at  a  given 
height  above  the  alfalfa.  One  can  conclude  that 
within  the  range  of  windspeeds  encountered, 
shearing  stress  is  almost  completely  dependent 
upon  the  physical  characteristics  of  the  surface. 


(cm./ctc.) 


0  100  200  300  400  500 

w160<IB.(cn»./*«c.) 

FuJi  Rt  19. — Wind  profile  parameters  for  alfalfa,  Julv 
29-31,  1960. 
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T \ble  1. — Functional  relationships  between  wimispeed  ami  roughness  length  (so),  ejfective  displacement 
(D),  ami  friction  relocity  (V*)  for  alfalfa  at  Ithaca ,  N.Y.,  July  29-31,  1960 


Profile  No. 


1 
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4  . 
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6  . 


7. 

8. 

9. 

10 

11 
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Date  and 

Height 

Velocity  | 

Zo 

D 

V*!k 

time 

E.t.l. 

Cm. 

Cm. /sec. 

Cm. 

Cm. 

Cm.  >sec. 

90 

87 

7/30/60 

1027 

100 

120 

107 
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15.  6 

-48.  4 

88.  8 

160 

175 

90 

118 

7/30/60 

1002 
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120 
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163 

25.  2 

-23.  3 

121.  2 
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205 

90 

130 

7/30/60 

0956 

100 

120 

145 

180 

28.  5 

-  18.  8 

140.  9 

160 

225 

7/30/60 

90 

1O0 

138 

156 

17.  8 

-32.  8 

118.  2 

1206 

129 

188 

160 

232 

7/29/60 

90 

140 

1840] 

1852 

100 

120 

160 

197 

-  16.  4 

-38.  8 

122.  8 

1904] 
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245 

7/29/60 

90 

158 

1916 

7/30/60 

100 

120 

182 

218 

11.  0 

-46.  4 

114.  7 

0950 

160 

268 

7/30/60 

1008 

1014 

90 

100 

120 

160 

175 

197 

230 

278 

11.3 

-39.  8 

117.  6 

90 

212 

7/29/60 

1846 

100 

120 

235 

275 

3.  5 

-57.  8 

94.  6 

160 

320 

90 

235 

7/29/60 

1605 

100 

120 

260 

295 

10.  7 

-33.  5 

141.  7 

160 

350 

90 

250 

7/31/60 

100 

280 

3.  2 

-59.  0 

110.  4 

1259 

120 

325 

160 

380 

90 

270 

7/31/60 

100 

305 

3.  2 

-57  9 

118.  1 

1311 

120 

347 

160 

408 

7/31/60 

90 

280 

12311 

100 

315 

.  5. 3 

-49.  4 

137.  6 

13531 

120 

355 

160 

420 

90 

305 

7/31/60 

100 

340 

2  H 

-57.  8 

125.  0 

1317 

120 

387 

160 

450 

| 

90 

312 

7/31/60  | 

'00 

350 

1  3 

tS5.  0 

196.  2 

1 247 

120 

395 

I 

H.0 

455 

7,31/60 

90 

345 

1219\i 

100 

3«5 

2.  7 

55.  8 

137.  1 

1253/ 

120 

430 

160 

500 
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When  one  observes  the  motion  of  the  alfalfa 
plants  in  the  field,  two  factors  readily  become 
apparent  As  the  wind  velocity  increases,  the 
leaves  orient  themselves  in  the  direction  of  the 
wind,  and  the  field  appears  to  consist  of  waves; 
i.e.,  “smoothing  out  effect.” 

Under  low  velocities,  the  plants  stand  nearly 
erect  while  each  leaf  nets  as  a  momentum  sink. 
The  momentum  loss  is  probably  neatest  in  the 
aren  of  the  upper  one-third  of  the  plants;  however, 
some  wind  movement  is  apparent  down  to  the 
ground  level.  As  wind  velocities  increase,  how¬ 
ever,  the  upper  leaves  orient  themselves  in  the 
direction  of  the  wind  flow  and  the  plants  overlap 
each  other.  Thus,  a  complete  surface  of  smooth 
leaves  is  approached  and  is  accompanied  bv  a  drop 
in  the  roughness  length,  z0.  For  the  effect  to 
occur,  the  plant  must  he  high  enough  to  provide 
a  complete  ground  cover  and  consist  of  enough 
foliage  to  create  the  “wave”  appearance  under 
higher  velocities. 

Wheat  Investigations 


average  height  of  the  wheat  was  approximately 
130  en:.  Most  of  the  growth  had  taken  place 
l)etween  the  upper  leaf  and  the  head.  Thus,  the 
center  of  mass  was  displaced  upward. 

Figure  21  and  t abb's  3  and  4  present  the  data 
showing  that  the  roughness  length,  Zo,  and  effec¬ 
tive  displacement,  D,  decreased  while  the  friction 
velocity,  U*,  remained  fairly  constant  with  in¬ 
creasing  windspeed  above  a  “critical”  level.  Evi¬ 
dently,  the  wheat  crop  on  July  22,  1960,  acted 
similarly  to  alfalfa.  The  reversal  in  trends  for  the 
wheat  crop  during  the  two  different  stages  of  de¬ 
velopment  has  to  lx*  attributed  both  to  a  change 
in  the  physical  characters  of  the  crop,  and  to  the 
diversity  of  range  of  windspeeds  encountered.  In 
the  earlier  stage  (table  2),  the  winds  were  slower 
and  the  crop  shorter,  thus  more  rigid.  This  would 
permit  some  waving  and  leaf  flutter,  but  little 
streamlining.  At  the  later  dates  under  the  higher 
wind  velocities  when  the  crop  was  taller,  stream¬ 
lining  probably  occurred  much  as  it  did  in  alfalfa. 
It  is  of  interest  to  point  out  that  the  trends  ob¬ 
served  in  the  w-heat  at  the  lower  windspeeds  on 


The  air  turbulence  studies  in  a  wheatfield  were 
made  on  an  open,  nearly  level  hilltop  (Mount 
Pleasant)  located  adjacent  and  to  the  north  of 
Ellis  Hollow  Experimental  site  described  else¬ 
where  (14).  The  field  has  the  approximate 
dimensions  of  200  by  200  meters.  As  in  the 
alfalfa  investigations,  measurements  were  made 
under  conditions  favoring  small  temperature 
gradients,  and  with  the  same  equipment. 

It  was  observed  earlier  tb  it  within  the  custom¬ 
ary  5-ininute  averaging  period  there  were  periods 
of  relatively  steady  wind.  In  an  effort  to  separate 
the  steady  flow  from  the  gusts,  a  series  of  30- 
second  sampling  periods  were  used.  When  a 
series  of  two  or  more  consecutive  30-second  periods 
rendered  similar  readings  from  the  cup  anemom¬ 
eters,  the  profile  was  considered  to  be  semi-steady. 
By  selecting  semi-steady  periods  the  semi-log 
plotting  of  the  wind  profiles  proved  to  be  much 
more  precise  than  when  tine  \Uc  averaging 
period  was  employed. 

On  the  first  day  of  study  (July  9,  1960),  the 
wheat  averaged  almost  100  cm.  in  height  at  a 
density  of  60  stnlks  per  square  foot.  The  heads 
were  10  cm.  long,  and  located  right  nbove  the 
upper  leaves.  Wind  profiles  were  taken  when  the 
wind  blew  from  a  northeasterly  direction.  As  can 
lie  seen  in  figure  20  and  table  2,  V*  and  za  in¬ 
creased  with  an  increase  in  wind  velocity.  Also, 
the  effective  displacement,  />,  decreased  when  the 
wind  velocity  nlx>ve  the  crop  increased. 

In  regard  to  these  measurements,  it  must  be 
mentioned  that  the  wind  velocities  did  not  attain 
High  enough  values  to  cause  appreciable  bending 
of  the  wheat  stalks.  A  different  picture  was 
obtained  July  22  and  23,  when  measurements 
were  taken  at  the  same  site  At  this  time,  the 
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Fun  rk  20. — Wind  profile  parameters  for  wheat,  July  9. 

19(50. 
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Table  2. — Wind  profile  data  jor  wheat  at  Ithaca, 
NY,  July  9,  I960 


1 

Profile 

No. 

Height 

Velocity  j 

V!k 

1) 

tr 

Cm. 

Cm. lire.  | 

Cm  nee. 

Cm. 

Cm. 

100 

180 

110 

225 

i . . 

130 

>  125.  3 

-82.  8 

170 

'  Is  'l 

l  250 

,  K731 

! 

^■7*1 

110 

345 

2 _ 

130 

>  108.  1 

-91.  9 

0.  788 

170 

485 

250 

320 

110 

390 

3 . 

130 

455 

>  195.  6 

170 

545 

250 

670 

110 

4  . 1 

130 

535 

}  266.  5 

-49.  4 

170 

250 

780 

(cm. /sac.) 
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Fioi're  21. —  Wind  profile  parameters  for  wheat,  July 
22-23,  1960. 


the  earlier  date  parallel  those  found  in  mature 
corn  at  higher  windspeed  in  our  1959  investiga¬ 
tions  (14).  A  few  values  for  com  of  V*,  z^,  and 
D  as  a  function  of  windspeed  are  given  later  in 
table  5.  They  are  in  agreement  with  the  1959 
trends. 


Table  3. — Wind  profile  data  jor  wheat  at  Ithaca, 
N.  Y.,  July  22,  I960 


Profile 

No. 

Height 

Velocity 

V*/k 

D 

Cm. 

Cm.ltcc. 

Cm.lsec. 

C  m. 

Cm. 

(  140 

380 

150 

410 

1 

{  170 

440 

210.  2 

-46.  6 

15.  0 

I  210 

500 

[  290 

585 

140 

440 

150 

470 

2 . . 

170 

500 

>  221.  6 

-58.  8 

11.  1 

210 

585 

290 

670 

140 

550 

150 

585 

3. . 

170 

640 

■  228.  7 

-82.  2 

5.  24 

210 

735 

290 

840 

140 

585 

|  150 

640 

4 _ 

<  170 

725 

■  285.  8 

-91.6 

6.  22 

1  210 

840 

l  290 

990 

(  140 

640 

150 

725 

5 _ 

170 

820 

213.  0 

- 1 18.  0 

1.09 

!  210 

950 

- 

{  290 

1080 

Figure  22  gives  a  qualitative  diagram  of  the 
functional  relationships  for  the  three  crops  for 
comparative  purposes.  The  trends  on  wind- 
speed  for  the  different  crops  are  chiefly  a  matter  of 
rigidity  of  momentum  sink  elements,  whereas 
the  order  of  magnitude  of  the  values  is  attributed 
to  the  momentum  sink  geometry.  The  alfalfa 
is  the  least  rigid,  the  taller  wheat  next,  followed 
by  the  shorter  wheat,  with  corn  being  the  most 
rigid.  A  maximum  Z0  probably  takes  place 
in  all  three  crops.  lx>w  enough  windspeeas  in 
alfalfa  were  not  encountered  to  clenrlv  define 
this  point,  whereas  it  might  take  very  high  winds 
to  observe  a  reversal  in  corn.  The  wheat  would 
fall  between  the  alfalfa  and  corn  in  characteristics. 
Of  course,  all  these  relationships  would  depend 
upon  stage  of  crop  maturity,  plant  density,  and  a 
host  of  other  factors  too  numerous  to  mention. 

TURBULENT  EXCHANGE  WITHIN  THE 
PLANT  CANOPY 

The  exchange  processes  within  the  vegetative 
canopy  are  theoretically  much  more  complex  than 
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Fini'RE  22. — Idealized  roughness  parameters  for  alfalfa,  wheat,  and  corn  as  a  function  of  windspeed. 


those  above,  because  of  the  complex  distribution  of 
sources  and  sinks  within  the  “crop  volume.” 
Instrument  and  interpretative  limitations  on  the 
experimental  level  are  restrictive,  too.  None¬ 
theless,  we  have  attempted  to  arrive  at  a  semi- 
quantitative  picture,  at  least,  of  some  transfer 
characteristics  within  the  wheat  and  corn  canopies. 

Referring:  again  to  equation  113,  we  have 


which  states  that  the  downward  flux,  Q ,  of  some 
entity  is  equal  to  t he  product  of  the  transfer  co¬ 
efficient,  K,  and  the  gradient  of  the  entity  across 
some  plane,  dc,dz.  If  we  assume  that  the  transfer 
coefficient  for  momentum  is  quantitatively  equal 
to  the  transfer  coefficient  for,  say  C02,  water 
vapor,  and  heat,  then  experimentally  we  need  to 
determine  the  transfer  coefficient  for  momentum 
from  windspeed  characteristics  and  the  gradient 
for  the  entity  in  question  (i.e.,  C02,  water  vapor 
or  heat),  in  order  to  calculate  the  exchange  rates 
for  the  latter. 

I  sing  Prandtl’s  model  for  rough  flow  in  a  pipe 
(33)  and  making  the  aSSuIl  iptions  that  turbulence 
is  fully  developed  and  isotropic  within  the  vege¬ 
tative  canopy,  we  have. 

K(z,)=ly!(w’)2  (121) 


where  K(zt)  is  the  momentum  transfer  coefficient 
at  nominal  height,  zt;  l  is  the  “mixing  length”; 
and  w'  is  the  vertical  windspeed  fluctuation,  or 
“vertical  eddy  velocity.”  The  eddy  velocity,  w' , 
is  defined  as  the  difference  between  the  instan¬ 
taneous  velocity,  w,  and  a  mean  velocity,  w. 

Experimentally,  w'  and  l  were  evaluated  from 
hot  wire  measurements  made  with  a  Hastings 
Model  HR-2  nondirectional  hot  wire  anemometer. 
The  Hastings  Model  HR-2  is  an  instrument  that 
measures  all  components  of  velocity  in  an  additive 
manner.  When  measurements  of  the  mean  hori¬ 
zontal  flow  were  made  at  a  point,  the  values  re¬ 
ceived  on  the  recording  chart  were  not  actually 
true  horizontal  values.  looking  at  the  flow 
pattern  in  two  dimensions,  the  hot  wire  measures 
tfie  mean  horizontal  flow  (vector  A)  and  the 
resolved  component  (vector  D)  of  the  horizontal 
fluctuation  (vector  B)  and  the  vertical  fluctuation 
(vector  C). 
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Table  4. — Wind  profile  data  jor  wheat  at  Ithaca , 
N.Y.,  July  23,  I960 


Profile 

Height 

Veloc- 

V*lk 

2o 

I) 

No. 

itv  1 

Cm. 

Cm.lsec. 

Cm.Jtec. 

Cm. 

Cm. 

25 

20.  7* 

50 

34  6* 

75 

41.  2* 

1. 

90 

70  f 

’  1 23  ±  5 

7.  2±.  9 

-77  ± 

too 

1 40 

120 

21 5 1 

160 

305 1 

240 

380  f 

25 

50 

37.  7* 

75 

48.  6* 

2 . 

90 

120t 

131-7 

5.  9  i  1 

—  75±2 

100 

195t 

120 

265 1 

160 

345 1 

240 

440 1 

25 

24  6* 

50 

52.  4* 

75 

54.  3* 

3 _ 

90 

140t 

160-6 

7.  4±.  8 

—  72±2 

100 

215f 

120 

295 1 

160 

400  f 

240 

500 1 

25 

25.  6* 

50 

51.  8* 

75 

64.  5* 

4 _ 

90 

1 80 1 

241  ±3 

14.  2±4 

-6C- 1 

100 

250  f 

120 

345 f 

160 

470  f 

240 

610f 

25 

29.  6* 

50 

89.9* 

75 

98.  0* 

5 . . 

90 

•205 1 
280  f 

100 

120 

370  f 

160 

490 1 

240 

640 1 

♦Hot  wire  data;  t  anemometer  cup  data. 


A  typical  hot  wire  trace  may  be  represented 
thus: 


The  mean  velocity  was  determined  by  averaging 
the  readings  at  intervals  of  two-thirds  second  over 
a  30-second  run.  If  the  horizontal  wind,  u,  is 


considered  in  a  steady  state  (a  to  b),  the  fluctu¬ 
ations  during  this  period  are  assumed  to  be  the 
fluctuations  of  the  vertical  wind  component,  w'. 
An  averaging  period,  however,  from  0  to  t  would 
include  a  large  eddy  regime.  Then,  averaging 
over  this  period  would  render  a  root  mean  square 
of  the  fluctuations  for  the  horizontal  component, 

■yfu’)*,  of  the  large  eddy.  On  the  other  hand,  by 
selecting  for  steady  periods,  a  good  estimate  of 

V( to ')*  may  be  obtained  (S3). 

Use  is  made  of  Taylor’s  Theorem  (34),  which 
employs  a  statistical  approach  to  turbulence. 
Using  an  entire  30-second  steadv  state  hot  wire 
record  and  displacing  it  two-thirds  of  a  second  for 
each  time  increment,  tu  we  assigned  the  whole 
record  equal  weight  in  an  autocorrelation,  Rt, 
where  Rt  is  the  correlation  between  w'h  at  the 
time  <i  and  u>J|+2/3  at  time  t1+2,3,etc.  A  good  es¬ 
timate  for  the  mixing  length,  l,  will  then  be: 


where  tn — time  for  Rt= 0. 

A  typical  autocorrelation  curve  is  shown  in 
figure  23.  It  is  noticed  that  instead  of  lim  Rt 

t-¥oo 

=0,  Rt  crosses  zero  and  then  oscillates  around  it. 
This  shows  that  the  ideal  case  is  approached  in 
this  evaluation,  but  never  achieved. 

Having  estimated  the  vertical  mixing  length, 
l,  and  the  vertical  component  of  velocity,  we  may 
then  determine  the  eddy  diffusion  coefficient  at 
any  point  in  space  from  equation  121.  Then. 
froni  this  relationship  and  the  wind  profile,  we 
may  estimate  the  shearing  stress  (momentum 
exchange)  by: 

t  —  pK(:()duldz,  or 

('(>2  exchange  P=K(zi)d  C02/ds,  (123) 

where  P  is  the  (’02  exchange  rate,  and  d  ('() 2-dz 
is  the  gradient  of  U02. 

Wheat  Investigations 

On  July  23,  t he  mast  was  placed  so  that  the  cup 
anemometers  were  at  heights  of  90,  100,  120,  160, 
and  240  cm.  above  t he  ground.  The  site  was  the 
same  ns  that  used  for  the  previous  measurements 
with  tlm  wheat  130  cm.  high  On  this  day,  how¬ 
ever,  a  Hastings  Model  HR  2  hot  wire  anemom¬ 
eter  with  an  0.5-second  time  constant  was  used  to 
obtain  mean  wind  velocities  at  heights  of  25,  50. 
and  75  cm.  in  the  plant  canopy. 

Since  there  was  only  one  hot-wire  instrument  to 
scan  the  chosen  individual  levels  in  the  vegetative 
canopy  of  t he  wheat,  a  “normalizing”  procedure 
had  to  be  adopted  to  permit  a  complete  mean 
windspeed  profile  evaluation.  This  consisted  of 
inking  numerous  semisteady  30-second  runs  with 
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FiGt  re  23. — Sample  autocorrelation  (Rt)  as_a  function  of 
laj;  time  (At)  times  mean  \vi!ids]>eed  (»<)  for  wheat, 
July  23.  1960. 

the  hot  wire  anemometer  at  a  given  level  in  the 
erop  simultaneous  to  numerous  semisteady  30- 
seeond  profile  runs  with  the  nip  anemometers 
above  the  erop.  Hot  wire  measurements  at  the 
various  levels  in  the  canopy  then  were  selected 
on  nearly  identical  “above  the  crop  profiles.” 
Five  complete  profiles  were  selected  in  this  way. 

Values  for  the  five  profiles  are  given  in  table  4 
and  figure  24.  There  are  three  observations  of 
immediate  importance:  (1),  the  velocities  at  25 
cm.  are  nearly  the  same  for  all  wind  regimes  cov¬ 
ered;  (2),  the  values  of  du/d:  between  50  and  75 
cm.  are  quite  small  for  each  profile;  and  (3),  the 
point  of  inflection  relative  to  the  crop  surface  ap¬ 
pears  to  be  within  20  cm.  below  the  crop  top. 

From  hot  wire  measurements  made  in  the  wheat 
(130  cm.  high)  and  cup  measurements  made  in  and 
above  the  wheat,  two  profiles  that  will  be  consid¬ 
ered  in  detail  in  this  section  are  plotted  in  figure 
25.  Using  the  method  just  described,  we  calcu¬ 
late  the  values  for  K(:t)  for  the  three  heights  of 
25,  50,  and  75  cm.  within  the  crop.  Above  the 
cn>p,  K(Zi)  was  calculated  from  t lie  semilog  profile, 
using  equation  120. 

Figure  26  shows  the  distribution  of  K(zt)  in¬ 
creasing  with  height  above  the  ground  and  m- 
creasicg  with  windspeed  above  the  crop.  The 
distribution,  K(z, ),  between  75  cm.  and  top  of  the 
crop  was  not  determined,  but  it  is  believed  to  be 


Z 

(cm.) 


V  (cm./i«c.) 

Fic.che  24. — Wind  profiles  in  wh^atfield,  July  23,  1960 

rather  eomplex,  based  upon  the  results  found  in 
corn. 

Two  items  of  importance  should  be  pointed  out. 
In  contrast  to  the  findings  reported  by  Penman 
and  Ixmg  (22)  for  wheat,  the  transport  coefficient 
here  is  considerably  larger  than  for  molecular  dif¬ 
fusion  (104  grenter)  nnd  the  wheat  crop  did  not 
appear  to  “seal  itself”  with  increasing  windspeed; 
i.e.,  the  transport  coefficient  within  the  crop  was 
sensitive  10  windspeed  above  the  erop. 

Calculations  for  shearing  stress  are  given  in 
figure  27.  In  this  case,  shear  does  appear  to  be 
considerably  attenuated  within  the  upper  half  of 
the  crop. 

Corn  Investigations 

On  August  2,  1960,  wind  profiles  were  taken  in 
Ellis  Hollow  cornfield  (14)  about  30  feet  north¬ 
west  of  the  central  instrument  line.  Northwest 
winds  blew  approximately  45°  to  the  corn  rows. 
The  sky  became  cloudy  after  2:30  that  afternoon, 
so  that  near  isothermal  conditions  existed  during 
the  test  period. 

Wind  profiles  were  taken  with  cup  anemometers 
above  tin*  corn  at  heights  of  250,  300,  and  350  cm. 
The  three  cup  anemometers  were  used  as  a  refer¬ 
ence  for  “normalizing”  the  hot  wire  nu  asurements 
of  25,  50,  75.  100,  150,  200,  230,  and  30)  cm.  above 
the  ground.  The  crop  averaged  alx.ut  240  cm. 
high  nnd  was  approximately  HO  percent  tasselcd. 
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1960. 


At  this  stage  of  development,  there  is  an  unusually 
high  concentration  of  leaves  at  the  top  of  the  plant. 
Subsequent  rapid  vertical  growth  of  the  stem  to 
full  height  eventually  leads  to  a  more  even  dis¬ 
tribution  of  leaves  along  the  stems. 

Four  wind  profiles  were  chosen  that  cover  the 
spectrum  of  winds  that  blew  that  afternoon.  Fig¬ 
ure  28  gives  the  graphical  analyses  of  the  turbulent 
boundary  layer  profiles  above  the  crop.  Too  few 
levels — three — were  taken  to  determine  adequately 
the  profile  parameters,  VT*,  z$,  and  D,  by  any 
statistical  method.  Nonetheless,  by  use  of  judg¬ 
ment  from  past  experience  in  the  graphical  method 
with  mature  com  (14)  realistic  values  were  arrived 
at  in  this  case.  Pertinent  data  are  given  in  table  5. 

The  measurements  taken  in  and  above  the 
crop  canopy  under  similar  winds  were  plotted 
in  figure  29.  The  solid  lines  were  drawn  in  as  a 
“best  fit”  approximation.  Three  points  should 
be  me.itioned  at  this  time: 

(1)  There  is  fairly  good  agreement  between 


Figure  26. — Transfer  coefficient  for  momentum  (Km)  as 
a  function  of  height  ( z )  in  and  above  a  wheat  crop, 
July  23,  1960. 


Table  5. — Wind  profile  data  taken  above  com  with 
cup  anemometers  at  Ithaca ,  N.Y.,  August  2 ,  1960 


Profile 

No. 

Height, 

r. 

Veloc- 
ivy,  u 

V*/k 

I) 

Cm. 

Cm. /sec. 

Cm. /sec. 

('m. 

Cm. 

250 

130 

1 

- 1 

300 

160 

56.  2 

r.  4 

-180 

1 

350 

ISO 

! 

250 

215 

) 

2..  i 

300 

3S5 

98.  5 

9.  6 

-  130 

350 

350 

J 

250 

300 

| 

3.. 

300 

405 

144.  8 

12  S 

-  100 

350 

435 

250 

420 

) 

4.. 

300 

470 

177.7 

17.  0 

-  80 

1 

350 

505 

the  cup  anemometer  and  hot  wire  at 
the  300-cm.  level; 

(2)  duldz  is  maximum  just  below  the  crop 
top; and 


z 

(cm.) 
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Figcre  27. — Drag  (r)  as  a  function  of  height  (*)  in  and 
above  a  wheat  crop,  July  23,  1960. 

(3)  du/dz  is  quite  small  between  25  and  100 
cm. 

In  explanation  of  the  last  observation,  it  was 
determined  that  the  vertical  component  of  velocity 
was  8  percent  of  the  mean  horizontal  component 
above  the  crop;  however,  this  value  increased 
as  the  height  of  the  measurement  was  decreased. 
At  the  25-crn.  position,  it  was  estimated  that  the 
vertical  component  was  about  30  percent  of  the 
mean  horizontal  flow.  Since  the  vertical  com¬ 
ponent  is  additive,  du/dz  between  25  and  100 
cm.  should  be  greater  than  shown  in  figure  29 
and  table  6.  Therefore,  the  previous  assumption 
of  isotropic  flow  is  less  correct  as  the  measure¬ 
ments  approach  the  ground. 

The  statistical  mixing  lengths  were  calculated 
at  each  point,  using  the  method  described  above. 
The  values  obtained  are  listed  in  table  7. 

The  values  for  above  the  crop  (300  cm.)  are  at 
best  approximate;  however,  they  are  in  the  right 
order  of  magnitude.  Since  the  ratio  of  l/z= 0.4 
within  the  turbulent  boundary  layer  (Karman), 


Z  ♦  0  (cm.) 


y  (cm./**c.) 

Figuhe  28. — Log  profiles  of  wind  above  cornfield,  August 
2,  1960. 

Table  6. — Wind  profile  data  taken  with  hot  wire 
anemometer  at  Ithaca ,  N.Y.,  July  2,  1960 


Height, 


Winds  peed, 


u,  according  to  profile  No. — 


(cm.) 


1 

2 

3 

4 

25-  _ 

Cm. /tec. 
12.  6 

Cm. /sec. 

18.  0 

Cm.lsec. 
21.  0 

Cm.  /tec. 
43.  0 

50... 

19  0 

20.  9 

30.  0 

48.  0 

75... 

19.  7 

26.  9 

32.  0 

.  _  .  ... 

100 

21.  5 

30.  0 

27.  4 

60.  1 

150  i 

23.  5 

44.  7 

52.  0 

100.  0 

•200... 

70.  4 

124.  0 

156.  7 

249.  5 

230 

111.  2 

166.  6 

228.  4 

300  0 

300 

158.  5 

266,  5 

425.  7 

488.  8 
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Table  7. — Mixing  length  l  (z){,  as  a  function  of 
height ,  at  Ithaca,  N.Y.,  Aug.  2,  1960 


Height, 
t>  (cm.) 

Mixing  length,  /,  for  profile  No. — 

1 

2 

3 

4 

Cm. 

Cm. 

Cm. 

Cm. 

25 _ 

21.0 

25.0 

25.  5 

30.  0 

60 _ 

31.0 

33.0 

36.0 

47.  0 

75 

39.  5 

45.  0 

47.  5 

100 _ 

50.5 

50.0 

53.0 

55.  5 

160 _ 

44.0 

53.  0 

78  0 

93.  0 

200 _ 

76.5 

75.  5 

100.0 

110.  0 

230 _ 

183.0 

196.0 

205.  0 

220.  0 

300 _ 

63.0 

74.  5 

140.  0 

170.  0 

one  would  expect  close  to  the  same  value  at  300 
cm.  Using  the  ground  as  the  reference  plane,  we 
obtain  l/z  range  from  0.21  to  0.565.  With  the 
zero-plane  displacement  given  as  the  reference, 
the  values  range  from  0.485  to  0.71. 

The  mixing  length  values  at  25  cm.  are  con¬ 
sidered  to  be  too  large.  This  is  attributed  to  the 
high  fraction  of  vertical  component  of  velocity 
when  calculating  the  mean  horizontal  flow. 
Within  this  region,  l/z  is  approximately  1. 

The  values  for  l/z  are  fairly  constant  between 
50  and  200  cm.  They  are  in  the  magnitude  of  0.5. 
However,  the  most  mysterious  values  were  ob¬ 
tained  at  a  point  10  cm.  below  the  crop  top 
(z=230).  At  this  point,  the  mixing  lengths  are 
larger  than  those  values  obtained  in  the  boundary 
layer  above  the  crop,  and  are  therefore  suspect. 
The  observation  made  while  the  data  were  being 
analyzed  indicates  that  the  reason  for  these  high 
values  originates  from  the  abnormally  long  period 
of  the  cycle  about  the  mean  veloc’'  \ 

At  the  position  of  230  cm  ,  the  frequency  of  one 
cycle  about  the  mean  was  from  three  to  five  times 
the  ones  previously  encountered.  It  is  suggested 
that  this  may  be  due  to  the  elastic  waving  of  the 
stalk  and  the  upper  leaves  superimposed  upon  the 
normal  frequency.  The  longer  cycles  were  respon¬ 
sible  for  a  much  less  rapid  drop  in  the  auto¬ 
correlation  with  time,  and  nence  a  larger  .alue  for: 


Rtdt. 


Values  for  the  eddy  diffusion  coefficient,  K(zt), 
were  calculated  from  equation  121  (table  K).  They 
were  then  plotted  in  figure  30.  The  dashed  lines 
in  figure  30  represent  K(:,)  as  calculated  from  the 
seindog  profile,  and  the  solid  lines  are  from  ‘he 
hot  wire  data.  The  shearing  stress  was  calculated 
from  equation  122  (table  9).  Figure  31  shows  the 
distribution  of  the  shearing  stress. 


Figure  29. — Wind  profiles  in  and  above  cornfield,  August 
2,  1960. 

Table  8. —  Transfer  coefficient,  K(zt),  as  a  function 
of  height  at  Ithaca.  N.Y.,  August  2,  i960 


Height, 

*.  (cm.) 

K(zi)  for  profile  No. — 

1 

! 

2 

3 

4 

Cm* /sec. 

Cm.*lsec. 

Cm.*  /sec. 

Cm.* /sec. 

25 _ 

94 

208 

298 

882 

50 

190 

287 

428 

1,349 

75. 

181 

486 

561 

100. 

359 

540 

1,039 

1,576 

150 

339 

1, 211 

2,  223 

2,  818 

200 

1, 306 

1, 910 

2,  700 

3,  894 

230 

5,  106 

5,  606 

6,  724 

7,  854 

300.  . 

1,  796 

2,712 

5,  250 

7,  208 

SUMMARY 

The  determination  of  Km  in  and  above  the 
momentum  sink  elements  using  hot  wire  anemom¬ 
eter  data  affords  an  independent  check  on  the 
semilog  profile  method  for  determining  Km,  and 
also  permits  a  distribution  analysis  of  where  the 
ip  Mnentuin  is  being  transferred  in  the  vegetative 
si _  Of  particular  interest  here  is  the  corn- 
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Z 


(cm.) 


Figvbe  30. — Distribution  of  transfer  coefficient  for  momentum  (Km)  in  and  above  cornfield,  August  2,  1960. 


Table  9. — Shearing  dress  as  a  function  of  height 
at  Ithaca ,  N.  Y. ,  August  2,  1960 


Height, 

Shearing  stress  for  profile  X 

o. — 

2,  (cm.) 

1 

> 

3 

4 

25 

Dynes/cm .* 
0.  038 

Dynes/cm.* 
0.  083 

Dynes/cm* 
0.  119 

Dynes/cm* 
0.  353 

50. 

.076 

.  115 

.  171 

.  540 

75. 

073 

.  194 

.  224 

100 

.  143 

.  227 

.  499 

.  946 

150 

203 

.  872 

2.  13 

3.  38 

200 

1.  88 

3.  21 

5.  51 

9.  35 

230  . 

9  19 

13.  45 

20.  17 

28.  27 

300 

1  08 

2.  28 

5.04 

8.  65 

parison  of  the  Km  values  derived  by  the  two 
methods  in  the  turbulent  boundary  layer  above 
the  corn  (z,  =  300  cm.)  in  figure  30.  In  three  out 
of  the  four  profiles  the  hot  wire  values  are  higher 
but  of  the  same  magnitude  as  the  log  profile 
values.  The  largest  divergence  occurs  at  the 
lowest  windspeed,  although  there  appears  to  be 
no  consistent  relation  between  divergence  and 
windspeed.  If  the  largest  divergence  at  slow 
windspeed  is  indeed  real,  it  should  be  anticipated 
from  the  following  factors: 

(1)  Rotating  cup  anemometers  are  subject 
to  greater  errors  at  slow  windspeeds ; 

(2)  There  are  errors  in  estimating  “free 
convection”  from  the  log  profile  caused 
by  thermal  gradients;  the  latter  are  more 
pronounced  at  low  windspeeds. 
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It  should  be  pointed  out,  however,  that  tem¬ 
perature  gradients  were  extremely  small  in  the 
turbulent  boundary  layer  over  the  com  (tempera¬ 
ture  gradients  may  have  been  significant  within 
the  canopy,  however),  and  that  tne  slowest  wind 
measured  above  the  crop  was  well  above  the 
stalling  speed  of  the  cup  anemometers. 

The  distribution  of  shear  in  figure  31  can  be 
considered.  Realizing  that  the  techniques  used 
and  the  assumptions  made  are  open  to  criticism, 
and  that  the  results  should  be  viewed  as  only 


Z  (cm.) 


semiquantitative,  we  should  point  out  that  the 
origin  of  shear  from  the  log  profiles  falls  on 
“planes”  between  80  and  180  cm.  above  the 
ground  (  —  l)  values.)  On  the  other  hand,  the 
hot  wire  evaluation  of  shear  within  the  canopy 
demonstrates  a  “volume”  distribution  of  shear 
in  a  complex  but  reasonable  fashion.  The  hot 
wire  quantitative  values  are  obviously  too  high 
if  the  assumption  is  made  that  the  total  shear  of 
the  “surface”  is  accurately  evaluated  by  the  log 
wind  profile  characteristics  above  the  “surface.” 


Figi  re  31. — Distribution  of  drag  (r)  in  and  above*  cornfield,  August  2,  19fi0. 


SUMMARY  AND  CONCLUSIONS 


1.  Two  theoretical  models  were  developed  for 
wind  flow  characteristics  in  plant  canopies  based 
upon  a  general  type  differential  equation  for  the 
wind  profile  both  inside  and  above  the  vegetative 
canopy  at  the  earth’s  surface. 

2.  A  testing  of  the  theoretical  models  depended 
upon  a  breakdown  of  the  wind  structure  into 
quasi-steady  and  transient  states.  It  was  assumed 
that  the  small  scale  eddies  caused  by  the  micro¬ 
structure  of  the  surfaces  (plant  surfaces)  are 
responsible  for  vertical  transport.  For  a  begin¬ 
ning,  ihe  simpler  case  was  studied  experimentally. 
Wind  flow  in  and  above  the  vegetation  at  the 
quasi-steady  state  produced  trends  predicted  by 
the  turbulent  flow  model. 

3.  The  turbulent  canopy  flow  model  postulated 
a  linear  growth  of  diffusivity  within  the  turbulent 
canopy  flow  layer  rather  than  the  log  velocity 
profile.  This  led  to  a  realistic  curvature  reversal 
of  the  velocity  profile  deep  in  the  canopy  flow 


layer.  This  must  be  a  basic  characteristic  of  the 
turbulent  wall  flow  within  the  roughness  elements. 

4.  Studies  of  the  log  velocity  wind  profile 
characteristic  over  various  kinds  of  vegetation 
revealed  a  complex  coupling  between  windspeed, 
roughness  length,  friction  velocity,  and  effective 
displacement  and  the  surface  geometric  and 
elastic  properties.  In  two  cases,  the  friction 
velocity  remained  constant  and  independent  of 
windspeed.  To  facilitate  these  and  future  studies, 
a  machine  computation  method  was  developed 
that  includes  an  evaluation  of  the  standard  errors 
of  the  logarithmic  wind  profile  parameters. 

5.  An  evaluation  of  the  transfer  coefficient  and 
shear  within  the  plant  canopy  revealed  a  reason¬ 
able  picture  of  their  spatial  distribution  within 
the  ‘‘crop  volume.”  This  serves  to  emphasize 
that  the  concept  of  the  roughness  length,  although 
very  useful,  is  in  reality  an  artifact. 
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